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Two  major  cost  components  in  distribution  system  are  inventory  and  transportation 
cost.  For  convenience,  the  inventory  and  transportation  decisions  are  determined  as 
two  disjointed  decisions  which  result  in  a high  system  total  cost.  Unifying  these  two 
decisions  can  lead  to  reduction  in  the  system  total  cost. 

This  dissertation  considers  the  problem  of  multiple  capacitated  warehouses  supplying 
multiple  retailers  with  stochastic  demand.  Its  purpose  is  to  develop  analytical  models 
and  solution  approaches  for  optimizing  a commodity  flow  in  a distribution  system  un- 
der uncertain  demand  based  on  the  sum  of  inventory  and  transportation  costs.  The 
research  aims  to  establish  improved  decision-making  methods  that  integrate  inventory 
replenishing  policy  and  transportation  strategy. 
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The  developed  model  formulation  extends  the  simple  single-period  stochastic  inven- 
tory problem  (newsboy  or  newsvendor  problem)  to  incorporate  the  transportation  cost 
into  the  objective  function.  In  addition,  the  capacity  constraints  are  added  to  impose 
over-supplying  at  warehouses.  Two  forms  of  the  transportation  cost  functions  are  con- 
sidered: linear  and  fixed  charge  cost  functions. 

The  Lagrange  multiplier  method  is  employed  to  solve  the  linear  transportation  cost 
model.  This  method  works  well  only  with  small  problems.  In  the  fixed  charge  cost 
model,  the  dynamic  slope  scaling  procedure  (DSSP)  is  employed  to  develop  the  scenario- 
based  heuristic.  The  DSSP  is  an  effective  and  efficient  approach  to  estimate  an  optimal 
solution  to  a fixed  charge  network  flow  problem.  Based  on  a set  of  test  problems,  com- 
putational results  show  that  the  DSSP  scenario-based  heuristic  obtains  the  solution  in 
reasonable  computational  time.  This  obtained  solution  to  the  approximated  problem 
(scenario-based  model)  is  used  to  evaluate  the  true  objective  value  in  the  original  non- 
linear problem,  which  gives  the  upper  bound  to  the  original  problem.  Since  the  optimal 
solution  to  the  original  problem  is  unknown,  the  conclusion  on  the  upper  bound  quality 
cannot  be  established.  Alternatively,  the  gap  between  upper  bound  and  lower  bound 
can  be  used  as  a criterion.  The  best  lower  bound  is  obtained  by  the  development  of 
Lagrangian  relaxation  model.  Computational  experiments  show  that  the  performance 
of  the  DSSP  heuristic  is  excellent,  based  on  a set  of  test  problems. 

Finally,  we  improve  the  DSSP  heuristic  by  employing  the  Lagrangian  relaxation 
method  to  solve  a linear  transportation  cost  and  eliminate  scenario  generating  work. 
Based  on  a set  of  test  problems,  the  Lagrangian  relaxation-based  heuristic  provides  a 
better  solution  in  less  computational  time  than  the  scenario-based  heuristic. 
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CHAPTER  1 
INTRODUCTION 

1.1  Background  and  Motivation 

Supply  chain  management  (SCM)  has  become  one  of  the  most  studied  research  ar- 
eas in  operations  research.  A supply  chain  involves  a series  of  activities:  purchasing, 
manufacturing,  stocking,  and  distributing.  This  activity  series  represents  a flow  of  com- 
modities from  procuring  raw  materials  to  delivering  to  end  customers.  Along  the  flow, 
one  or  more  facilities  of  the  following  (vary  from  company  to  company)  are  involved: 
factories,  distribution  centers,  warehouses,  and  retailers.  Thus,  SCM  is  the  management 
of  this  flow  which  is  equivalent  to  logistics  management. 

Embedding  in  distributing  activities,  inventories  and  transportation  are  two  key  el- 
ements in  the  logistics.  Since  the  future  is  not  certain,  inventories  are  served  as  precau- 
tionary assets.  They  protect  companies  against  commodity  shortage  which  gradually  de- 
stroys customer  satisfaction,  company’s  reputation,  and  company’s  market  position.  As 
for  transportation,  customers  are  not  typically  located  in  the  same  places  as  manufactur- 
ers. Thus,  transportation  is  a crucial  element  in  the  supply  chain  to  move  commodities 
to  the  hand  of  customers  which  is  the  closing  point  of  the  commodity  flow. 

The  combined  inventory  and  transportation  cost  is  equal  to  one-half  to  two-thirds 
of  total  logistics  costs  [8].  Robert  V.  Delaney  reported  in  “State  of  Logistics  Report” 
that  in  1999  American  business  logistics  costs  were  $921  billion  (9.9  percent  of  nominal 
Gross  Domestic  Product).  This  figure  had  increased  about  7 percent  from  1997.  Since 
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the  inventory-transportation  cost  plays  an  important  role  in  the  logistics  total  cost, 
there  is  an  increasing  need  for  an  efficient  and  effective  distribution  strategy  in  which 
the  company’s  logistics  total  cost  can  be  reduced.  Determining  such  a strategy  requires 
a well-developed  analytical  model  and  analysis  techniques  that  coordinate  inventory 
management  and  transport  strategy. 

The  purpose  of  this  dissertation  is  to  develop  analytical  models  and  solution  ap- 
proaches for  determining  an  efficient  distribution  strategy  in  a system  of  multiple  ware- 
houses supplying  multiple  retailers.  The  developed  model  aims  to  optimize  a commodity 
flow  so  that  the  system  total  cost  is  minimized.  The  scope  encompasses  transportation 
cost  of  shipping  commodities  from  warehouses  to  retailers  via  different  transport  chan- 
nels (connection  between  a warehouse  and  a retailer),  and  inventory  decisions  (holding 
and  shortage)  at  retailers.  In  this  introductory  chapter,  a brief  overview  of  a single 
period  stochastic  inventory  problem  and  a transportation  problem  is  presented  in  Sec- 
tion 1.2  and  1.3.  The  organization  of  the  remainder  of  this  dissertation  is  framed  in 
Section  1.4. 


1.2  The  Single  Period  Stochastic  Inventory  Problem 
As  its  name  suggested,  this  problem  deals  with  the  one-period,  periodic  review  in- 
ventory problem  with  a probabilistic  demand.  A classical  example  of  this  type  is  the 
newsboy  or  newsvender  problem.  It  has  the  simplest  model  structure  which  considers 
only  a single  commodity,  no  fixed  ordering  cost,  zero  initial  inventory,  and  no  constraint. 
In  this  problem,  a newsboy  must  order  newspapers  to  sell  at  the  beginning  of  a day.  If 
the  newsboy  carries  too  many  newspaper,  then  many  worthless  newspapers  are  left  at 
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the  end  of  the  day.  On  the  other  hand,  a newsboy  who  orders  too  few  newspapers  will 
lose  opportunity  to  make  additional  profit.  The  question  is  how  many  newspaper  should 
be  ordered  in  order  to  maximize  the  newsboy’s  profit. 

Applications  of  this  static  one-shot  type  decision  often  appear  in  the  manufacturing 
of  highly  seasonal  products  such  as  toys,  high-fashion  clothes,  and  sports  event  souvenirs. 
Similar  to  the  newsboy  problem,  the  products  may  become  obsolete  at  the  end  of  the 
season  if  the  production  is  exceeding  the  demand.  On  the  other  hand,  a loss  of  profit 
may  be  induced  by  shortage  of  products  resulting  from  a limited  production.  The 
management  must  wonder  what  would  be  the  optimal  production  level  so  as  to  minimize 
their  expected  total  cost  or  maximize  their  expected  profit.  An  important  key  of  this 
decision-making  process  is  the  uncertainty  or  randomness  of  demand.  If  the  demand  is 
known,  then  this  decision  would  become  easier  to  make. 

The  newsboy  problem  and  its  solution  are  presented  in  many  introductory  textbooks 
in  operations  research  [54,  59],  operations  management,  or  inventory  management  [52, 
47].  The  closed- form  formulations  for  computing  the  expected  total  cost  are  presented 
in  Lau  [36].  These  formulations  are  derived  with  various  demand  distributions.  The 
inventory  replenishing  policy  for  the  newsboy  problem  is  known  as  an  order-up-to  policy 
resulting  from  the  tradeoff  between  holding  (overstocking)  and  shortage  ( understocking) 
cost.  Usually,  the  optimal  order-up-to  level  is  denoted  by  S*.  This  problem  is  more 
complicated  when  a fixed  cost  is  incurred  for  replenishment  but  not  otherwise.  This 
fixed  cost  may  include  an  order  processing  cost,  an  administrative  cost,  and  so  on. 
Taking  into  account  this  fixed  cost,  the  replenishing  policy  is  modified  to  (s,S)  policy. 
Under  this  policy,  when  the  initial  inventory  level  at  the  beginning  of  the  period  is  less 
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Expected  Total  Cost 


Figure  1.1:  Expected  cost  function  for  newsboy  model  with  fixed  ordering  cost 

than  s,  the  order  is  placed  to  bring  its  level  to  S.  The  {s,S)  policy  incorporates  two 
decision  variables:  s and  S.  The  variable  s specifies  when  to  make  an  order,  while  both 
variables  s and  S specify  how  many  units  to  order. 

Let  G{»)  and  K denote  the  expected  cost  function  and  fixed  ordering  cost  respec- 
tively. The  optimal  s*  is  derived  from  the  equation  below. 

G(s*)  = K + G{S*) 

where  S*  can  be  derived  by  the  same  method  in  the  order-up-to  policy.  Figure  1.1  shows 
that  it  is  better  not  to  order  when  the  inventory  level  falls  between  (s*,  S*).  Otherwise, 
the  fixed  ordering  cost  would  increase  the  total  expected  cost.  As  K decreases,  s* 
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becomes  larger.  At  the  point  where  K is  zero,  s*  is  equal  to  S*.  It  is  noted  that  the 
order-up-to  policy  is  a special  case  of  (s,S)  policy  in  which  s is  equal  to  S. 

Many  studies  in  the  literature  have  been  done  to  extend  the  simple  newsboy  problem. 
The  investigation  of  multi-commodity  or  multi-product  type  system  appears  in  several 
references  [54,  45,  37,  38].  Some  studies  focus  on  incorporating  constraints  (resource 
restrictions)  [52,  45,  37,  38].  Many  contributions  are  devoted  to  the  demand  ‘distribution 
free’  approach.  This  approach  uses  the  information  of  estimated  mean  and  variance 
without  the  knowledge  of  demand  distribution  [18,  26,  44,  43].  Its  goal  attempts  to 
decrease  the  value  of  missing  information  (demand  distribution). 

1.3  The  Transportation  Problem 

The  transportation  problem  (TP)  is  a special  case  of  the  minimum  cost  network 
flow  problem  (MCNFP).  The  MCNFP  determines  the  least  cost  of  distributing  the 
commodities  through  a network  such  that  demands  at  certain  nodes  are  satisfied  by 
available  supplies  of  other  nodes  [2].  Different  from  the  basic  MCNFP,  the  TP  has  a 
bipartite  network  structure.  That  means  the  node  set  is  partitioned  into  two  subsets: 
demand  node  subset  and  supply  node  subset.  A supply  node  is  connected  to  a demand 
node  by  an  arc.  The  flow  on  each  arc  has  only  one  direction  from  the  supply  node  to  the 
demand  node.  Usually  the  demands  are  known  and  the  supplies  are  finite  at  the  supply 
node.  One  of  applications  of  the  TP  is  the  distribution  of  products  from  warehouses 
to  retailers.  Suppose  that  a system  consists  of  m warehouses  supplying  n retailers.  Let 
i index  the  warehouses  in  the  set  / = {!,...  ,m}  and  j index  the  retailers  in  the  set 
J = {1,...  ,n}.  A flow  on  each  arc  incurs  a unit  variable  cost,  Cjj.  Let  Xij  be  the  amount 
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shipped  from  warehouse  i to  retailer  j.  The  available  capacity  at  the  warehouse  i is  Si, 
while  the  amount  required  at  the  retailer  j is  dj.  Then  the  mathematical  formulation  of 
the  TP  can  be  given  as  follows: 


n m 


mi”  EE 

j=l  i~l 

Cij  3Cij 

n 

s.t.  Xij 

3-^ 

lA 

< 

(Supply  constraints) 

m 

2=1 

> Sj  Vj 

(Demand  constraints) 

Xij 

> 

o 

Al 

(Nonnegativity  constraints) 

The  objective  is  to  minimize  total  cost  of  distributing  commodities  from  supply  nodes 
to  demand  nodes  subject  to  all  constraints.  In  supply  constraints,  the  total  amount 
shipped  from  each  warehouse  cannot  exceed  its  capacity.  In  demand  constraints,  the 
total  amount  received  at  each  retailer  must  satisfy  its  demand.  The  last  constraint 
ensures  that  all  flows  cannot  be  negative.  This  problem  can  be  solved  efficiently  (see 
[2]  for  several  algorithms  for  transportation  problems). 

Another  form  of  the  TP  is  a fixed  charge  transportation  problem  (FCTP).  Now  the 
cost  flow  in  each  arc  has  two  components:  fixed  and  variable.  The  fixed  component  does 
not  depend  on  the  amount  shipped,  but  rather  on  whether  or  not  there  is  a shipment 
flow  on  that  arc.  In  the  logistics  context,  the  fixed  component  may  include  the  minimum 
freight  charge  for  delivery  of  the  commodities  or  an  order  processing  charge.  The  FCTP 
belongs  to  the  class  of  the  minimum  concave  cost  network  flow  problem  (MCCNFP) 
which  is  an  A/”P-hard  problem.  This  means  there  is  no  known  algorithm  which  can  solve 
this  class  of  problem  in  polynomial  time.  Guisewite  and  Pardalos  [30]  survey  general 


characteristics  of  the  MCCNFP  and  its  solution  methods. 
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Transformed  into  a zero-one  mixed  integer  programming  problem,  the  FCTP  can  be 
solved  using  branch  and  bound  algorithms  to  obtain  an  exact  optimal  solution.  Many 
studies  have  been  done  to  enhance  the  existing  branch  and  bound  based  algorithms.  One 
recent  work  by  Bell  et  al.  [11],  develops  a capacity  improvement  technique  to  improve 
their  branch  and  bound  algorithm. 

1.4  Organization  of  the  Dissertation 

The  material  presented  in  this  dissertation  is  organized  into  six  chapters. 

• Chapter  2 deals  with  the  problem  of  deriving  the  optimal  inventory  replenishing 
policy.  It  begins  with  a review  of  various  stochastic  inventory  problems  of  both 
single-echelon  and  multi-echelon  systems.  Then,  it  develops  a single-period  two- 
echelon  stochastic  model  with  a fixed  ordering  cost  for  a single-warehouse  single- 
retailer system.  The  comprehensive  analytical  result  is  presented. 

• Chapter  3 introduces  transportation  cost  in  a capacitated  multi-warehouse  multi- 
retailer single-period  stochastic  inventory  problem.  A review  of  the  literature  in 
minimizing  inventory-transportation  cost  model  is  presented  at  the  beginning.  A 
nonlinear  programming  model  is  developed  to  minimize  the  inventory  and  trans- 
portation costs  simultaneously.  This  model  facilitates  in  determining  an  integrated 
inventory  replenishing  policy  and  a transporting  plan.  Developed  first,  the  linear 
transportation  cost  model  applies  the  Lagrange  multiplier  method  for  analyzing 
the  optimal  integrated  policy. 

• Chapter  4 substitutes  the  fixed  charge  transportation  cost  in  the  optimization 
model  presented  in  Chapter  3.  The  fixed  charge  model  optimal  value  is  estimated 
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by  the  DSSP  heuristic  through  a scenario-based  approximated  problem.  The  com- 
putational results  for  a wide  variety  of  problem  sizes  are  reported. 

• Chapter  5 measures  the  quality  of  the  upper  bound  obtained  from  Chapter  4 by 
comparing  it  with  a lower  bound.  The  best  lower  bound  can  be  obtained  from  the 
development  of  a Lagrangian  relaxation  model.  This  relaxation  model  employs  the 
subgradient  optimization  algorithm  to  compute  the  lower  bound.  The  extensive 
comparison  results  are  presented. 

• Chapter  6 develops  a new  Lagrangian  based  DSSP  heuristic.  This  approach  em- 
ploys the  Lagrangian  relaxation  method  to  solve  the  linear  transportation  cost 
problem.  As  a result,  it  eliminates  an  approximation  of  demand  distribution 
through  the  scenario  generating  process. 

• Chapter  7 describes  directions  for  future  research.  It  surfaces  a broad  variety  of 
research  challenges  to  be  treated  in  much  greater  depth. 


CHAPTER  2 

THE  SINGLE-PERIOD  TWO-ECHELON  STOCHASTIC  INVENTORY  PROBLEM 

2.1  Introduction 

Most  inventory  systems  in  a snpply  chain  are  of  a multi-echelon  type,  the  man- 
agement of  which  is  complex  and  challenging.  With  this  type  of  system,  dependency 
among  echelons  cannot  be  ignored.  In  practice  individual  stages  often  minimize  cost 
based  on  local  information  to  develop  simple  models  and  policies.  This  approach  may 
not  minimize  the  total  cost  of  the  whole  system.  Research  has  also  shown  that  policies 
based  on  local  information  can  lead  to  distorted  and  amplified  demand  requirements  at 
higher  echelons  [39].  Thus  the  dependence  among  echelons  should  be  considered  as  a 
coordinated  system,  which  we  would  like  to  operate  at  minimum  cost. 

While  many  works  have  been  devoted  to  the  multi-echelon  inventory  control  problem, 
systems  with  uncertain  demands  have  received  attentions  only  in  recent  years.  Thus, 
this  chapter  embarks  on  the  multi-echelon  stochastic  inventory  problem  with  analyzing 
a simple  two-echelon  inventory  problem  with  probabilistic  demand.  The  simple  newsboy 
problem  is  modified  to  incorporate  an  additional  echelon  inventory  and  fixed  ordering 
cost  is  present. 

The  material  of  this  chapter  is  organized  as  follows.  Section  2.2  provides  a review 
of  literature  on  supply  chain  management.  Section  2.3  describes  the  framework  of  the 
single-period  two-echelon  stochastic  inventory  problem.  Section  2.4  introduces  notation, 
presents  the  mathematical  formula,  and  derives  the  optimal  policy.  Finally,  Section  2.5 
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illustrates  the  derivation  of  the  optimal  policy  under  the  negative  exponential  demand 
distribution. 


2.2  Literature  Review 

Recently,  many  researchers  have  turned  their  attention  to  supply  chain  management 
(SCM)  [56].  A multi-echelon  inventory  system  is  one  of  the  most  interesting  problem 
in  SCM.  In  fact,  the  multi-echelon  inventory  system  has  been  studied  since  the  1960s. 
The  most  classical  work  is  that  of  Clark  and  Scarf  [21].  Their  model  focuses  on  a 
serial  multi-installation  inventory  system,  and  an  optimal  policy  is  determined  for  serial 
systems  using  a recursive  computation  approach.  While  the  demand  is  originated  at  the 
lowest  installation,  no  setup  cost  exists  in  any  installation  except  the  highest  installation. 
The  term  “echelon”  stock  is  introduced.  The  echelon  stock  is  referred  to  the  stock  at 
any  given  installation  plus  stock  in  transit  to  or  on  hand  at  downstream  installations. 

The  probabilistic  multistage  serial  inventory  system  is  extended  by  Chen  and  Zhen  [20] 
to  include  setup  costs.  For  echelon-stock  (r,  nQ)  policies,  they  propose  an  efficient 
heuristic  for  determining  near-optimal  control  parameters  in  multistage  serial  systems 
with  compound  Poisson  demands.  Their  algorithm  starts  with  determining  both  lower 
and  upper  bound  on  the  cost  function  by  over-  and  under-charging  a penalty  cost  to 
each  upstream  stage  for  holding  inadequate  stock.  Then  these  bounds  are  minimized 
to  obtain  a heuristic  solution,  that  is  near  optimal.  In  addition,  they  also  provide  an 
algorithm  that  ensures  the  optimal  solution.  However,  this  algorithm  requires  more 
computational  effort  than  is  practical  using  existing  computing  power. 
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Love  [40]  applies  dynamic  programming  recursion  to  develop  an  algorithm  for  finding 
an  optimal  nested  schedule,  for  an  n-period  single-product  inventory  model.  The  model 
is  A^-facility  serial  inventory  system  with  deterministic  demand  and  separable  concave 
production  and  storage  costs.  The  demand  occurs  at  the  lowest  echelon  (labeled  N). 
The  nested  schedule  implies  that  if  there  is  a production  at  facility  j for  a given  period  i, 
then  all  its  successors  must  produce  in  period  i.  To  guarantee  the  existence  of  a “nested” 
optimal  schedule,  the  storage  costs  are  not  decreasing  as  the  product  flows  down  stream 
and  the  production  costs  are  not  increasing  in  time. 

A relatively  sophisticated  model  of  a two-level  inventory  system  of  a single-item 
source  supplying  n parallel  warehouses  with  probabilistic  demand  and  no  setup  cost  is 
investigated  by  Gross  [29].  This  model  does  not  require  the  same  demand  distribution 
for  all  warehouses.  Transshipping  between  warehouses  is  considered,  because  it  might  be 
more  economical  to  transfer  inventory  from  overloaded  warehouses  to  understocked  ones 
instead  of  ordering  from  the  central  source.  The  model  neglects  the  lead-time  between 
placing  and  receiving  an  order  and  that  for  initiating  and  receiving  transshipment.  Thus 
the  total  expected  system  costs  are  associated  with  variable  ordering  cost,  transshipment 
cost,  and  holding  and  shortage  costs.  To  minimize  the  expected  total  system  costs,  an 
iterative  procedure  is  applied  to  a one-period  model  to  provide  a practical  solution. 

The  two-level  inventory  system  with  a single  source  supplying  n parallel  depots  is 
extended  by  Mahoney  and  Sivazlian  [41]  to  incorporate  setup  costs  and  interdependence 
among  the  depots.  Under  a joint  ordering  policy  or  (a,  S)  policy,  the  demand  at  each 
depot  is  assumed  to  have  an  Erlang  distribution,  and  no  transshipment  between  depots 
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is  allowed.  Given  an  order-up-to  level,  Si's,  and  the  “reorder”  level,  an  expression  for 
steady  state  stockout  probability  at  each  depot  is  derived. 

Ernst  and  Pyke  [22]  examine  a two-level  distribution  system  composed  of  a warehouse 
and  a retailer,  which  faces  random  demand.  The  model  considers  both  inventory  and 
transportation  costs.  Thus  the  expected  total  cost  function  is  defined  as  holding  and 
stockout  costs  at  the  warehouse  and  the  retailer,  and  the  transportation  cost  for  shipping 
goods  from  warehouse  to  retailer.  An  algorithm  is  given  that  determines  the  optimal 
base  stock  policy  at  the  retailer  and  the  warehouse,  the  optimal  in-house  and  contracted 
truck  capacity,  and  the  optimal  review  period  length. 

The  study  of  a three-level  distribution  system  appears  in  Chan  and  Simchi-Levi  [19]. 
Their  system  consists  of  a single  outside  vendor,  a fixed  number  of  warehouses  and  many 
geographically  dispersed  retailers.  The  study  offers  an  efficient  algorithm  for  three  level 
distribution  systems  so  that  long-run  average  cost  is  as  small  as  possible.  The  long-run 
average  cost  is  defined  as  the  total  inventory  holding  and  transportation  costs  (from 
the  outside  vendor  to  the  warehouses  and  from  the  warehouses  to  the  retailers).  The 
distribution  strategy  does  not  keep  stock  at  the  warehouses,  a strategy  often  referred  as 
“cross-docking.”  In  other  words,  as  soon  as  a fully  loaded  truck  from  the  vendor  arrives 
at  each  warehouse,  the  warehouse  loads  the  goods  onto  other  trucks,  which  deliver 
consolidated  shipments  of  various  items  to  retailers.  Each  retailer  faces  a deterministic 
demand  rate,  which  is  assumed  identical  over  all  retailers.  No  shortage  or  backlogging 
is  allowed.  Initially,  each  retailer  is  assigned  to  a single  warehouse.  This  assignment  is 
similar  to  a bin-packing  problem  where  the  capacity  of  the  warehouse  is  analogous  to 
the  bin  capacity  and  retailers  are  equivalent  to  items  that  require  packing  in  the  bins. 
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Retailers  assigned  to  the  same  warehouse  require  the  same  delivery  schedule.  Retailers 
at  each  warehouse  are  then  partitioned  into  clusters.  Retailers  in  the  same  cluster  receive 
shipments  on  the  same  truck.  The  partitioning  problem  is  similar  to  finding  an  optimal 
traveling  salesman  tour. 

According  to  Lee,  Padmanabhan,  and  Whang  [39],  a demand  variation  phenomenon 
under  stochastic  environment  in  a serial  supply  chain  is  called  “bullwhip”  effect.  In 
general,  small  variations  or  nearly  steady  consumer  demand  are  amplified  into  high 
variation  at  higher  echelon  in  the  chain.  By  investigating  the  amplification  of  inter- 
echelon order  variation  in  multi-echelon  supply  systems.  Burns  and  Sivazlian  [16]  suggest 
two  sources  that  create  this  amplification:  (1)  a legitimate  and  unavoidable  inventory 
adjustment  and  (2)  an  unwarranted  “false  order”  effect.  A new  decision  rule  is  proposed 
to  suppress  false  orders  and  avoid  stockouts. 

As  reviewed  in  Sivazlian  [53]  and  Sivazlian  and  Stanfel  [54],  literature  relevant  to  our 
analysis  offers  optimal  ordering  strategies  and  some  characteristics  of  a multi-commodity 
inventory  systems  with  probabilistic  or  uncertain  demands  over  a single  period.  Interde- 
pendence is  assumed  among  the  demands  of  commodities  as  well  as  the  set-up  costs.  In 
the  analysis,  the  marginal  distributions  of  each  commodity  are  only  needed  rather  than 
the  joint  density  distribution.  The  analysis  is  initiated  with  a two-commodity  inventory 
system.  The  authors  then  extend  their  analysis  to  an  m-commodity  inventory  system 
(m  = 1,  2,...).  In  an  m-commodity  system,  the  number  of  possible  ordering  decisions 
equals  2”*,  with  an  order  being  placed  in  2"^  — 1 of  these  cases.  Thus  the  objective  is  to 
determine  the  minimum  expected  total  cost  among  those  ordering  decisions  in  a single 
period. 
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Figure  2.1:  A single-commodity  two-echelon  system  in  series 
2.3  Problem  Description 

Consider  a single-period  single-commodity  inventory  system  consisting  of  three  levels 
operating  in  series  as  shown  in  Figure  2.1.  The  system  is  composed  of  one  retailer, 
one  warehouse,  and  one  source  with  infinite  supplies.  During  a period,  depletion  of 
inventory  at  the  retailer  is  induced  by  customer  demand  (D),  while  inventory  depletion 
at  the  warehouse  {Dw)  is  only  induced  by  requests  from  the  retailer.  Lead-time  at  both 
warehouses  is  neglected.  For  simplicity,  we  assume  that  the  demand  at  warehouse  equals 
demand  at  retailer  {Dw  = D).  Let  (poi*)  be  a density  function  of  customer  demand 
and  let  P>d{u)  du  and  D = up>d{u)  du. 

We  also  assume  that  the  customer  demand  must  always  be  met.  In  other  words,  the 
unfulfilled  demand  is  backordered,  at  some  penalty  cost  per  unit  associated  with  it.  One 
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important  assumption  to  be  made  is  to  use  echelon  stock  level  instead  of  installation 
stock  level.  The  echelon  stock  at  the  warehouse  is  defined  as  the  stock  at  installation 
warehouse  plus  the  stock  at  echelon  retailer.  The  echelon  stock  at  the  retailer  is  the 
same  as  installation  stock  because  it  is  at  the  lowest  echelon  of  the  system. 

2.4  Model  Formulation  and  Analytical  Results 
Now  we  define  some  notations. 

Notation 

cr  variable  procurement  cost  at  the  retailer  in  $/(unit) 
cw  variable  procurement  cost  at  the  warehouse  in  $/(unit) 

Kr  fixed  order  cost  at  the  retailer  in  $ 

Kw  fixed  order  cost  at  the  warehouse  in  $ 

h*ji  holding  cost  at  installation  retailer  in  $/[(unit)(unit  time)] 

Hr  holding  cost  at  echelon  retailer  in  $/ [(unit)  (unit  time)] 

holding  cost  at  installation  warehouse  in  $/[(unit)(unit  time)] 
h\v  holding  cost  at  echelon  warehouse  in  $/[(unit)(unit  time)] 
shortage  cost  at  installation  retailer  in  $/[(unit)(unit  time)] 

Pr  shortage  cost  at  echelon  retailer  in  $/[(unit)(unit  time)] 

pIy  shortage  cost  at  installation  warehouse  in  $/[(unit)(unit  time)] 

Pw  shortage  cost  at  echelon  warehouse  in  $/[(unit)(unit  time)] 
xr  inventory  level  prior  to  making  a decision  at  echelon  retailer 
{xr  can  be  negative) 
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Xw  inventory  level  prior  to  making  a decision  at  echelon  warehouse 
(x\v  can  be  negative) 

zr  quantity  received  at  echelon  retailer  {zr  > 0) 
zw  quantity  received  at  echelon  warehouse  {zw  > 0) 

Ur  inventory  level  following  a decision  at  echelon  retailer 

{VR  = Xr  + Zr) 

yw  inventory  level  following  a decision  at  echelon  warehouse 
{Uw  — + Z\y) 

For  this  particular  system,  when  a decision  for  replenishing  stock  is  to  be  made,  there 
are  4 possible  courses  of  action  or  ordering  decisions  available.  These  are 

1.  order  for  the  retailer  only, 

2.  order  for  the  warehouse  only, 

3.  order  for  both  warehouses,  and 

4.  do  not  order  anything. 

2.4.1  Cost  Structure  of  the  Problem 

We  now  develop  a cost  function  for  inventory  that  includes  procurement  cost  and 
holding  and  shortage  costs. 

Procurement  cost 

For  each  warehouse,  the  two  basic  components  of  the  procurement  cost  are  a fixed 
cost  and  a variable  cost.  The  fixed  cost  is  incurred  only  if  an  order  for  replenishment 
is  placed.  The  variable  cost  depends  on  the  quantity  received  at  each  warehouse.  Let 
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P{xji,xw)  denote  the  procurement  cost  function;  then 

P{xr,  Xw)  = k{zR,  Zw)  + CrZr  + CwZw  (2.1) 

where  k{zR,  Zw)  represent  the  fixed  portion  of  the  procurement  cost,  and  crZr  and  cwZw 
denote  the  variable  procurement  costs  at  the  retailer  and  the  warehouse  respectively.  The 
fixed  portion  of  the  procurement  cost  may  include  the  delivery  cost  such  as  a minimum 
freight  charge,  the  driver’s  wages,  the  cost  of  processing  an  order,  and  so  forth.  The 
variable  procurement  cost  is  the  cost  per  unit  of  commodity  and  delivery,  excluding  the 
fixed  portion.  We  assume  that  the  variable  cost  per  unit  at  the  warehouse  is  strictly  less 
expensive  than  at  the  retailer.  That  is  cr  > cw-  In  general  k{zR,zw)  can  take  many 
values  depending  on  the  aforementioned  ordering  decisions.  Thus  we  can  define  it  as 
follows: 


k{zR,  z\v)  — KrSr  + Kw^w  (2.2) 

where  5r  = 1 ii  zr  > 0,  otherwise  0,  and  = 1 if  zvy  > 0,  otherwise  0. 

Holding  and  shortage  costs 

At  the  end  of  a period,  there  may  be  an  amount  of  stock  left  over  or  a shortage 
of  inventory  at  either  or  both  warehouses.  The  total  expected  holding  and  shortage 
cost  incurred  over  the  period  at  both  warehouses  is  a function  of  the  echelon  stock  level 
(yR,yw)  following  the  decision  process  at  the  start  of  the  period:  let  L{yR,yw)  denote 
the  expected  holding  and  shortage  costs  at  the  warehouse  and  the  retailer.  Thus  the 
magnitude  of  L{yR,  yw)  is  directly  effected  by  the  initial  echelon  inventory  level  {xr,  Xw) 
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and  the  ordered  amount  {zr^zw)-,  which  varies  based  on  the  ordering  decision.  Since  we 
use  echelon  stock  as  our  decisions,  the  holding  cost  and  shortage  cost  should  be  charged 
based  on  echelon  costs.  The  echelon  holding  and  shortage  costs  at  the  warehouse  are 
defined  as  follows: 


h\Y  — h^r 
Pw  =P*W 

Since  a part  of  holding  and  shortage  cost  at  the  retailer  have  already  been  charged  at 
the  warehouse.  Thus,  the  echelon  holding  and  shortage  costs  at  the  retailer  are  defined 
as  follows: 


Pr  = P*r-P*w 


Now,  we  can  express  L{yR,yw)  as  follows. 


rVR  poo 

L{yR,yw)  = hR  iVR-  OvniOd^  + PR  / i^-yR)^D{Od^ 
Jo  J yn 


-\-h\v 


rVw  roo 

/ {yw  - + Pw  / i^-yw)pD{Od^  (2-3) 

Jo  J yw 


The  function  L{yR,yw)  is  separable  in  j/r  and  yw-  This  means  that  we  can  express 
L{yR,  yw)  as  the  sum  of  a function  of  yR  only  and  of  a function  of  yw  only: 


L{yR-,  yw)  — L{yR)  + L{yw) 


(2.4) 
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where  L{yji)  and  L{yw)  can  be  expressed  as  follows: 


L{yR)  = hR  {vr  - + PR  (2.5) 


L{yw)  — hw  {Vw  - 0pd{0  + Pw  I ~ yw)pD{C)  d^  (2.6) 


It  is  interesting  to  note  that  both  L{yR)  and  L{yw)  are  twice  differentiable  at  all  points 


yR  and  yw  > ^ for  which  v?d(«)  is  continuous.  We  note  that  L"{yR)  = {Hr  + PR)(pD{yR.) 
and  L''{yw)  = {hw  + Pw)PD(yw)-  Thus,  L{yR,yw)  is  strictly  convex  in  yR  and  yiy. 

2.4.2  Objective  Function 

Now  we  wish  to  determine  an  optimal  replenishment  policy  so  as  to  minimize  the 
sum  total  of  the  procurement  cost  P{xr,xw)  and  the  expected  holding  and  shortage 
L{yRiyw)  costs  during  a single  period.  Let  C{xr,xw)  denote  the  minimum  total  ex- 
pected cost  for  the  period  when  following  an  optimal  policy.  Then  C{xr,  xw)  can  be 
expressed  as  follows. 


C{xr,xw)  — 


mm 

yR>XR,yw>xw 


{P{xr,  Xw)  + L{yR,  yw)} 


yR>XR,yw>xw 


min  {k{yR  - xr,  yw  - xw)  + CR{yR  - xr) 


-\-cw{yw  — Xw)  + L{yR,  yw)}- 


(2.7) 


Define 


G{yR,  yw)  — CRyR  cwyw  + L{yR,  yw)- 


(2.8) 
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Then 

C{xr,xw)^  {Kvr- ^R^Vw  - xw)  - crXr- cwxw + G{yR,yw)]-  (2.9) 


The  function  G{yR,yw)  plays  an  important  role  in  determining  optimal  ordering 
decisions.  As  defined  by  relation  (2.8),  it  is  evident  that  if  L{yR,yw)  is  strictly  convex, 
so  will  be  the  function  G{yR,yw)  ■ As  defined,  G{yR,yw)  will  possess  an  absolute 
minimum  occurring  at  a point  {Sr,Sw),  such  that  at  yR  = Sr  and  yw  = Sw  where 
dG{yR,  yw)/yR  = 0 = dG{yR,  yw)/yw-  The  existence  of  {Sr,  Sw)  is  guaranteed  because 
of  the  strictly  convexity  of  G{yR,  yw)- 

Now  we  want  to  find  Sr  and  sw-  From  (3.6)  and  (2.8),  we  can  express  G{yR,  yw)  as 
follows: 


yR>^R,yw>xw 


G{yR,  yw)  — CRyR  + cwyw  + L{yR)  + L{yw) 


or 


G{yR,  yw)  — G{yR)  + G{yw)- 


(2.10) 


Once  we  know  {Sr,  Sw),  we  can  find  Sr  and  Sw  from  the  equations  below. 


G{sR,yw)  — Kr-\- G{SR,yw) 


G{yR,sw)  — Kw G{yR,Sw) 


or 


G(s«)-G(5r)  = Kr 


(2.11) 


G{sw)  — G{Sw)  — Kw 


(2.12) 
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Figure  2.2:  Possible  regions  for  values  of  inventory  level  prior  to  making  decision  at  the 
retailer  and  the  warehonse  (xw) 

The  uniqueness  of  sr  and  sw  are  guaranteed  because  G{yR)  and  G{yw)  are  strictly 
convex  function  of  yji  and  yw 

To  get  more  insight  into  this  problem,  we  can  divide  the  possible  solutions  of  xr  and 
Xw  into  4 regions  according  to  predefined  ordering  rules.  These  regions  are  represented 
by  Figure  2.2.  One  more  assumption  needs  to  be  made  at  this  point:  the  warehouse 
places  and  receives  an  order  before  the  retailer  does. 

• If  {xr^xw)  falls  in  area  I [xr  < sr,  and  sw  < xw  < •S'w))  the  retailer  will  place 
an  order  of  Sr  — Xr  in  order  to  bring  the  inventory  level  up  to  Sr.  Unfortunately, 
sometimes  the  warehouse  may  not  have  enough  inventories  to  fill  the  request. 
Thus,  there  are  two  likely  events  occurring  under  this  area: 
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- If  xw  > (Sr  — xr)  ; the  quantity  received  is  Sr  — xr  . We  will  call  this  action 
ordering  decision  1(a).  The  total  expected  cost  associated  with  this  decision 
can  be  expressed  as 

Cl[a){xR,Xw)  = Kr  — CrXr  — CwXw  + G{Sr,X\v)  (2-13) 

- If  Xw  < {Sr  — xr),  then  the  retailer  will  have  xw  after  ordering  rather  than 
Sr.  We  will  call  this  action  ordering  decision  1(b).  The  total  expected  cost 
associated  with  this  decision  can  be  expressed  as 

Cl{b){xR,  Xw)  = Kr  — CrXr  — CwXw  + G{xw,  Xw)  (2-14) 

• If  {xr,  Xw)  falls  in  area  II  {xr  < sr,  and  xw  < svy),  the  retailer  and  the  warehouse 
will  place  an  order  of  Sr  — Xr  and  Sw  — xw  in  order  to  bring  the  inventory  level 
up  to  Sr  and  Sw  respectively.  Thus,  the  amount  to  be  received  at  the  retailer 
is  Sr  — Xr,  and  at  the  warehouse  is  Sw  — xw-  We  will  call  this  action  ordering 
decision  2.  The  total  expected  cost  associated  with  this  decision  can  be  expressed 
as 


G2{xr,  Xw)  — Kr  + Kw  — CrXr  — CwXw  + G{Sr,  Sw)  (2.15) 

• If  {xr,xw)  falls  in  area  III  {sr  < xr  < Sr,  and  xw  < sw))  the  warehouse 
will  place  an  order  of  Sw  — Xw  in  order  to  bring  the  inventory  level  up  to  Sw- 
As  mentioned  above,  the  source  has  infinite  supplies,  so  the  warehouse  receives 
inventories  in  an  amount  of  Sw  — xw  ■ We  will  call  this  action  ordering  decision 
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3.  The  total  expected  cost  associated  with  this  decision  can  be  expressed  as 

Csixji,  xw)  = K\v  — crXr  — cwxw  + G{xr,  ^vi^)  (2.16) 

• Finally,  if  {xr,xw)  falls  in  area  IV  {sr  < xr,  and  sw  < xw),  both  warehouses 
will  not  order  anything.  We  will  call  this  action  as  ordering  decision  4.  The  total 
expected  cost  associated  with  this  decision  can  be  expressed  as 

Ci{xR,  Xw)  = —crXr  — cwXw  + G{xr,  Sw)  (2-17) 


2.5  An  Example  with  Negative  Exponential  Demand  Distribution 
Suppose  that  the  demand (D)  for  the  commodity  at  the  retailer  has  negative  expo- 
nential distribution  with  mean  Let 

ipviO  = (2.18) 

$o(^)  = 1 - A > 0,  0 < ,^  < oo.  (2.19) 


From  equation  (2.5),  (2.6),  (2.8),  and  (2.10),  we  can  write  G{yR)  and  G{yw)  as  follows: 
G{yR)  = {Hr  + pr)  {yR  - + y - (ph  - cr)^/?  (2.20) 

pyw 

G{yw)  = {hw  + Pw)  J {Uw  — O'Poi^)  d^  + — — {pw  — cw)yw-  (2.21) 

Then, 

— = 0 = {Hr  -h  PR.)^D{yR)  — {PR  — Cr) 


dG{yw) 

dyw 


— 0 — {hw  + Pw)^D{yw)  — {pw  — Cw) 
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Hence,  if  pr  > cr  and  pw  > cw,  Sr  > 0 and  Sw  > 0 can  be  obtained  from 


^d{Sr) 


Pr  - Cr 
hR  -\-pR 


(2.22) 


Pw  — Cw 
hw  + Pw 


(2.23) 


Substituting  (2.19)  in  (2.22)  and  (2.23),  we  obtain  Sr  and  Sw  as  follows: 


1 Hr+Pr 

Sr  = yIii-t — ^ 

A nR  + Cr 


Sw  — T In 


1 , hw  + Pw 


A hw  d"  (^w 


Now,  for  Sr  and  sw,  substitute  (2.18)  in  (2.20)  and  (2.21): 

f>~^yR  2 T)n 

G{yR)  = {hR  +PR){yR  H X ~ ~ 


(2.24) 


g-Ayw  I 

G{yw)  = {hw  +Pw){yw  H A^  ~A — ~ cw)yw-  (2.25) 


Substituting  (2.24)  in  (2.11)  and  (2.25)  in  (2.12),  we  will  obtain  Sr  and  sw  by  solving 
the  following  equations. 


g— g-AS/{ 

Kr  = {hR  +Pr)[{sr  H — ) — {Sr  H — )]  — {pr  — Cr){sr  — Sr) 


g— As(v  g— ASw 

Kw  = {hw  -hPw)[{sw  H — ) — {Sw  H — )]  — {pw  ~ cw){sw  — Sw) 


CHAPTER  3 

THE  LINEAR  TRANSPORTATION  COST  MODEL 
3.1  Introduction 

Two  major  cost  components  in  logistics  are  inventory  and  transportation  costs.  In 
general,  the  products  or  commodities  are  manufactured  at  the  factory  and  kept  in  the 
warehouse  before  being  shipped  out  to  retailers  or  end  customers.  The  transportation 
cost  is  incurred  by  moving  commodities  over  space  while  the  inventory  cost  may  be 
incurred  from  holding  excess  commodities  at  the  retailers.  The  transportation  cost  may 
be  reduced  by  shipping  commodities  in  larger  shipment  with  less  frequencies.  However, 
the  inventory  holding  cost  may  increase.  Determining  inventory  and  transportation  de- 
cisions separately  generally  lead  to  a high  system  total  cost.  As  a result,  there  is  an 
increasing  need  of  integrated  models  which  coordinate  inventory  control  and  transporta- 
tion planning  aspects  [12].  The  objective  is  to  determine  an  integrating  decision  so  as 
to  minimize  the  system  inventory  and  transportation  costs.  The  integrating  decision 
results  from  the  tradeoff  between  inventory  cost  and  transportation  costs. 

This  chapter  develops  a mathematical  model  to  minimize  inventory  and  transporta- 
tion cost  in  multi-warehouse  multi-retailer  system  with  stochastic  demand  subject  to 
finite  capacity  at  warehouses.  Initially,  we  investigate  the  model  of  linear  transporta- 
tion cost.  Since  the  objective  function  is  the  sum  of  the  inventory  convex  cost  function 


25 


26 


and  the  linear  transportation  cost,  the  problem  becomes  a nonlinear  optimization  prob- 
lem. Silver  and  et  al.  [52]  derive  the  solution  of  the  similar  problem  with  one  resource 
through  the  use  of  Lagrange  multipliers. 

The  chapter  is  organized  as  follows.  A literature  reviewed  is  presented  in  Section  3.2. 
Section  3.3  gives  details  of  the  problem,  introduces  necessary  assumptions  and  notation, 
and  presents  the  mathematical  formulation.  The  Lagrange  multiplier  method  is  de- 
scribed in  Section  3.4.  Then  the  analysis  of  the  two-warehouse  two-retailer,  single  ware- 
house multi-retailer,  and  multi-warehouse  multi-retailer  system  is  shown  in  Section  3.5, 
3.6,  and  3.7  respectively.  Finally,  we  conclude  the  chapter  in  Section  3.8. 

3.2  Literature  Review 

The  purpose  of  this  section  is  to  survey  the  broad  variety  of  integrated  inventory  and 
transportation  models.  As  reviewed  in  the  literature,  the  transportation  cost  is  defined 
in  different  circumstances.  It  may  refer  to  a function  of  cost  per  trip  [55,  17],  cost  per 
shipping  unit  (embedded  in  inventory  model),  or  vehicle  routing  cost  [23].  The  inventory 
cost  is  defined  based  on  deterministic  or  stochastic  demand  models.  In  the  deterministic 
model,  the  inventory  cost  refers  to  the  holding  cost  charged  against  excess  inventory  in 
the  system.  In  the  stochastic  demand  models,  in  addition  to  the  holding  cost,  there  is 
a shortage  or  penalty  cost.  It  is  charged  for  not  being  able  to  complete  the  demand. 

Early  studies  in  an  integrated  inventory-transportation  system  is  done  by  Federgruen 
and  Zipkin  [23].  They  develop  a stochastic  model  which  combine  the  inventory  allocation 
problem  and  the  standard  vehicle  routing  problem.  An  inventory  allocation  problem 
determines  optimal  inventory  replenishing  policies  among  several  locations  (e.g.,  retailers 
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or  end  customers)  by  sharing  a limited  resource  (e.g.,  finite  capacity  warehouse  or  fieet 
vehicles) . The  vehicle  routing  problem  ( VRP)  determines  a delivery  plan  (vehicle  routes) 
for  a fleet  of  vehicles  at  minimum  system  total  cost.  The  authors  consider  a single- 
period system  of  one  central  depot  multiple  locations  with  stochastic  demand.  They 
manipulate  the  model  in  order  to  partially  decompose  it  into  an  inventory  problem  and 
several  VRP  subproblems.  The  interchange  heuristic,  a well  known  solution  method  for 
the  deterministic  VRP,  is  modified  to  manage  their  stochastic  model. 

Recently,  Bertazzi  and  Speranza  [12]  have  surveyed  the  deterministic  models  and 
algorithms  for  the  minimization  of  inventory  and  transportation  during  the  last  two 
decades.  The  models  and  algorithms  are  classified  in  continuous  and  discrete  models. 
The  reviewed  models  are  considered  under  the  following  five  basic  types  of  logistics 
networks:  the  single  link  case  (simplest  network),  the  sequence  of  link  case,  the  one 
origin-multiple  destinations  case,  the  multiple  origins-one  destination  case,  and  the  mul- 
tiple origins-multiple  destinations  case.  Most  of  literature  in  the  continuous  models  are 
‘succinct  modeling’.  The  results  in  many  of  these  models  are  on  the  basis  of  the  eco- 
nomic ordering  quantity  {EOQ).  One  disadvantage  of  the  continuous  time  model  is  the 
solution  may  not  be  realistic  from  the  practical  point  of  view.  The  aim  of  Anily  and 
Federgruen  [3]  is  to  determine  the  inventory  replenishing  policy  which  is  incorporated 
with  the  transportation  plan  over  the  infinite  horizon.  They  develop  a single-item,  de- 
terministic model  in  one  origin-multiple  destinations  network.  The  depot  (an  origin)  is 
a transshipment  point  (i.e.,  there  is  no  inventory  kept  at  the  depot).  They  propose  a 
heuristic  exploiting  the  regional  partitioning  strategy  to  cluster  the  destinations  (retail- 
ers) into  regions.  It  is  possible  that  the  same  retailer  belongs  to  several  regions.  That 
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means  a retailer  can  be  assigned  to  different  routes  simultaneously.  Once  the  regions  are 
established,  an  efficient  route  within  each  region  is  determined.  The  derived  lower  bound 
and  upper  bound  are  compared  to  show  that  these  bounds  become  asymptotically  tight 
when  the  number  of  destinations  approaches  infinity.  Later  in  Anily  and  Federgruen  [4] 
this  model  is  extended  to  coordinate  the  inventory  at  the  depot  level.  The  inventory 
cost  is  computed  by  a more  complicated  function. 

The  single  link  case  (one  origin  one  destination)  discrete  time  deterministic  model 
is  considered  in  Speranza  and  Ukovich  [55].  The  objective  is  to  determine  the  trip  fre- 
quencies in  a multi-product  system  with  a given  set  of  feasible  frequencies.  That  means 
only  one  frequency  is  assigned  to  each  product.  In  their  model,  several  products  have 
to  be  shipped  by  vehicles  of  given  capacity  on  a common  link  from  an  original  node  to  a 
destination  node.  Four  different  situations  are  analyzed  based  on  the  shipping  policies 
and  consolidating  rules.  There  are  two  schemes  in  shipping  policies  that  each  product 
can  follow.  It  can  be  either  shipped  in  only  one  frequency  or  partially  shipped  with 
multiple  frequencies.  The  way  that  products  are  loaded  on  the  vehicles  is  determined  by 
the  consolidating  rules.  It  can  be  either  frequency  or  time  consolidation.  Each  situation 
is  formulated  in  a mixed  integer  linear  programming  model.  They  have  proven  that  the 
model  with  multiple  frequencies  and  time  consolidation  is  equivalent  to  the  model  with 
multiple  frequencies  and  frequency  consolidation.  The  numerical  example  shows  that 
the  best  scheme  is  allowing  prodncts  to  be  partially  shipped  (multiple  frequencies)  and 
share  the  same  truck  with  the  same  frequency  product. 

Burns  et  al.  [17]  use  an  analytical  approach  to  minimize  the  inventory  and  trans- 
portation costs  in  one  origin-multiple  destinations  network.  Their  model  considers  only 
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one  product  with  deterministic  demand.  Two  distribution  strategies  are  developed,  eval- 
uated and  compared.  The  first  strategy  is  direct  shipping  (separate  loads  directly  to 
each  customer).  The  second  strategy  is  peddling  (more  than  one  customer  in  one  load). 
In  the  peddling  strategy,  customers  are  divided  into  delivery  regions.  Each  truckload 
is  delivered  to  customers  within  the  same  region.  The  transportation  cost  for  the  ped- 
dling involves  three  stages:  line-haul  (from  the  supplier  to  the  first  customer  stop),  local 
delivery  (from  the  first  customer  stop  to  the  last  customer  stop),  and  back-haul  (from 
last  customer  stop  to  the  supplier).  While  direct  shipping  uses  the  EO(5-based  optimal 
shipment  size,  the  peddling  shipment  size  is  a full  truck.  In  addition,  the  results  denote 
that  the  peddling  is  superior  to  the  direct  shipping  with  an  increase  of  the  following  ele- 
ments: item  value,  distance  from  supplier,  customer  density  in  the  region,  and  inventory 
holding  cost. 


3.3  Problem  Descriptions  and  Model  Formulation 
This  section  describes  the  capacitated  inventory  and  transportation  problem  in  a 
multi-warehouse  multi-retailer  system  with  stochastic  demand. 

First  let  us  revisit  the  single-period,  single-commodity  uncapacitated  inventory  prob- 
lem with  stochastic  demand.  The  decision  to  be  made  is  how  much  to  order  at  the  be- 
ginning of  the  period  while  the  expected  total  cost  is  minimizes.  The  expected  total  cost 
is  comprised  of  two  components:  holding  cost  and  penalty  or  shortage  cost.  The  holding 
cost  accrues  as  a result  of  having  positive  inventory  while  the  penalty  cost  or  shortage 
cost  is  charged  against  a negative  level  of  inventory  (unable  to  satisfy  demands)  at  the 
end  of  period.  Since  the  objective  function  is  convex  (shown  in  previous  chapter)  and 
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there  are  unlimited  resources,  then  there  is  a unique  optimal  solution.  Let  x*  denote 
the  optimal  solution.  Then  a;*is  derived  from 


P 

h + p 


where  p is  a cost  per  unit  of  unsatisfied  demand  and  h is  a cost  per  unit  of  having 
positive  inventory  at  the  end  of  the  period. 

As  mentioned  in  Chapter  1,  in  a transportation  problem,  a set  of  supply  nodes 
serves  a set  of  demand  nodes.  The  decision  to  be  made  is  how  many  units  to  be  shipped 
on  each  arc  between  supply  nodes  and  demand  nodes  in  order  to  minimize  the  total 
transportation  cost.  This  decision  must  not  exceed  the  limited  resources  available  at 
the  supply  side  and  must  satisfy  the  amount  needed  at  the  demand  side. 

Now  consider  a system  consisting  of  m warehouses  supplying  n retailers  (see  Fig- 
ure 3.1).  Let  i index  the  warehouses  in  the  set  I = {1,...  ,m}  and  let  j index  the 
retailers  in  the  set  J = {1, . . . ,n}.  Random  demands  take  place  at  retailers.  Ware- 
houses have  finite  capacity.  The  shipping  flow  can  be  induced  from  any  warehouse  to 
any  retailer.  When  the  inventory  model  is  extended  to  embrace  transportation  cost,  the 
optimal  solution  is  not  only  influenced  by  the  holding  and  penalty  cost  but  also  the  unit 
transportation  cost. 


3.3.1  Assumptions  and  Notation 

This  subsection  states  the  necessary  assumptions  and  notation. 
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Figure  3.1:  A capacity  constrained  multi- warehouse  multi-retailer  system 
Assumptions 

1.  The  capacitated  inventory  and  transportation  problem  is  formulated  as  a discrete- 
time, single  period,  multi-warehouse,  multi-retailer  problem. 

2.  The  warehouses  hold  single- type  commodities. 

3.  There  is  no  transshipment  between  warehouses  or  retailers. 

4.  The  shipping  decision  is  made  at  the  beginning  of  the  period. 

5.  The  decision  specifies  the  shipping  amount  and  its  associated  shipping  arc. 

6.  Supply  at  warehouses  is  limited. 

7.  Demands  at  retailers  are  stochastic  with  known  means  and  known  variances,  and 
independently  distributed.  Unsatisfied  demands  are  completely  lost  sales. 

8.  At  the  beginning  of  the  period,  the  inventory  level  at  all  retailers  is  zero. 
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9.  The  total  cost  includes  inventory  cost  (holding/penalty  cost)  and  transportation 
cost. 

10.  Inventory  costs  occur  at  retailers. 

11.  The  transportation  cost  can  have  a linear  form  (considered  in  this  chapter)  and  a 
fixed  charge  form  (considered  in  the  subsequent  chapters). 

Notation 

The  following  notation  is  used  in  the  rest  of  this  dissertation. 

I a set  of  warehouses 

J a set  of  retailers 

m number  of  warehouses 

n number  of  retailers 

Wi  warehouse  i in  set  I 

Rj  retailer  j in  set  J 

Xij  shipping  amount  (decision  variable)  from  Wi  to  Rj 

<PDj{Cj)  demand  probability  density  function  at  Rj 

cumulative  demand  distribution  function  at  Rj 
where  $d,(0)  = fo'  ipDj{u)du 
Dj  mean  of  demand  at  Rj 

Pj  penalty  or  shortage  cost  at  Rj 

hj  holding  cost  at  Rj 

Ci 


capacity  at  Wi 
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3.3.2  Model  Formulation 
The  involved  costs 

As  mentioned  earlier,  the  expected  total  cost  includes  the  inventory  cost  and  the 
transportation  cost.  Let  ICj  denote  the  inventory  cost  at  Rj,  and  Tij  denote  the  trans- 
portation cost  from  Wi  to  Rj.  It  should  be  noted  that  the  inventory  cost  is  evaluated 
from  the  the  total  amount  received  at  each  retailer  which  is  supplied  from  at  least  one 


rable  for  each  arc  ij  (see  Equation  (3.1)).  The  expected  inventory  cost  for  each  retailer 
is  expressed  as  follows. 


While  the  inventory  cost  is  a convex  function,  the  transportation  cost  can  take  any  form, 
eg.  linear,  fixed  charge,  etc.  Later  in  this  chapter,  we  consider  the  simplest  case,  a linear 
transportation  cost.  Now  let  TC{xij)  denote  the  expected  total  cost  of  the  system.  Then 
we  can  write  TC{xij)  as  follows. 


m 


warehouse.  Thus  the  inventory  cost  is  a function  of  ^^Xij.  That  means  it  is  not  sepa- 

i=0 


m 


(3.2) 
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Substitute  Equation  (3.1)  in  Equation  (3.2)  we  obtain 


TC{xij)  = ~ ) ¥^Dj  {^j ) 


j=0 


i=l 


poo  rn  \ 

+Pj  / ^ (^j  ~ ^^^ij)^Dj{^j)d^j  I 


+ Tij{Xij) 

j=0  i=0 


Model  formulation 

The  objective  is  to  find  the  optimal  integrating  inventory  and  transportation  decision 
that  minimizes  the  expected  total  cost. 

n / m 

min  hj  (J2  ^ij  - Q^Pdj  {Qd(j 

j=l  V i=l 

m \ 

+Pj  ^ij  ~ ^ ^ij)'fDj  i^j)d^j  j 

n m 

+ ^^Tij{xij)  (3.3) 

j=l  i=l 
n 

s.t.  "^^Xij  < Ci  Vi 

t=i 

Xij  >0  Vi,j 

There  are  two  sets  of  constraints.  The  first  set  is  the  capacity  constraints.  The  total 
amount  shipped  from  Wi  to  any  retailers  must  be  at  most  the  available  capacity  at  W,. 
The  second  set  is  nonnegativity  constraints. 

A linear  transportation  cost  model  formulation 


In  this  chapter,  we  consider  a linear  transportation  cost.  Let  Vij  denote  the  cost 
associated  with  moving  one  unit  of  commodity  from  the  Wi  to  the  Rj.  Then  Tij(xij)  is 
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expressed  as 


Tij(Xij^  — (3-4) 

Substitute  Equation  (3.4)  in  the  problem  (3.3),  the  formulation  becomes 

n / m 

hj  (Y  - QVd,  {Qd^j 

V i=l 

m \ 

+Pj  Xij  “ Y1  {^j)d^j  I 

t=l  / 

n m 

j=l  i=l 

n 

s.t.  Xij  < Ci  Mi 

J=1 

Xij  > 0 Mi,j 

3.4  Lagrange  Multiplier  Method  for  Inequality  Constraints  Problem 
This  section  describes  the  Lagrange  multiplier  method  to  solve  a nonlinear  program- 
ming problem  with  inequality  constraints.  The  method  is  summarized  from  [13,  14]. 
Consider  a problem  involving  inequality  constraints  (ICP) 

min  f{x) 

s.t.  gi  < 0,  . . . , Qrn  < 0 

Assume  that  / : 5R"  — >■  5?,  5'  : — > 5?”*  are  continuously  differentiable  functions.  Let  us 

introduce  additional  variables  cti, . . . , a^.  We  can  convert  the  ICP  above  to  the  equality 
constraint  problem  as  follows: 

min  f{x) 

s.t.  gi+aj  = 0,  . . . , gm  + al,  = 0 
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Let  X*  be  a local  minimum  for  our  original  problem  (ICP).  Then  (x*,q;*)  is  a local 
minimum  for  converted  problem,  where  a*  = {a\, , o;^), 

for  i = 1, . . . , m. 

For  all  j’s,  then  there  exist  unique  Lagrange  multipliers  such  that 

m 

V/(i')  + ^A,V!7,-(x-)=0, 

j = l 

Qj  + q;|  = 0 
2Xjaj  = 0 

3.5  The  Two- Warehouse  Two-Retailer  Problem 
This  section  analyzes  a simple  system  consisting  of  two  warehouses  and  two  retailers. 
We  can  formulate  the  problem  as  follows: 

min  hi  (pD^{^i)d^i 

+Pi  /i7i+x2i  (^1  ~ (^11  + ^12))  PDii(i)dCi 

+Ciia:ii  + C21X21 

+h2  ((xi2  + X22)  - 6)  PD2{i2)d^2 

+P2  I^2+X22  (^2  - (^12  + 2:22))  PD2i^2)d^2 

-\-V12X12  + V22X22 

s.t.  X11+X12  < Cl 

X21  + X22  ^ C2 


Xn,Xl2,X2l,X22  > 0 
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Converting  the  two  inequality  constraints  to  the  equality  case,  we  obtain 

^11  + ^12  + Oix  = Cl 

3^21  + 3;22  + ~ C2 

The  Lagrangian  function  of  the  problem  can  be  expressed  as 
F(a:ii,  a;i2,  ^21, 3:22,  Ai,  A2,  ai,  0.2)  = 

^Ill+X21  poo 

h\  / ((^11  + a;2i)  - +Pi  / (6  ~ (^11  + a:i2))<^Di(6)c?6 

■'O  '^Xll+X21 

pXl2+X22 

-\-ViiXii  + ^213^21  + ^2  / ((^12  + ^22)  — ^2)^D2i^2)d^2 

Jo 

POO 

+P2  / (^2  — {^12  + ^22))PD2{.^2)d^2  + ^123;i2  + '^22^22 

J XI2+XOO 


+Ai(a;ii  + a:i2  — C\)  + A2(3;2i  + X22  + <^2  ~ ^2) 


(3.6) 


By  rearranging  terms  and  using  integration  by  parts  technique,  Equation  (3.6)  is 
equivalent  to 

F(a;n,  a:i2, 0:21,  X22,  Ai,  A2,  Q!i,  0:2)  = 

r>Xii+X2l 


^DA^i)d^i  +PiD\  - [pi  - Vii)xii  - (pi  - U2i)a;2i 

'•X12+X22 


r 

(hi+Pi) 

Jo 

+ {^2  + P2)  / ^D2  (6)c?6  + P2D2  - (P2  - ^^12)2:12  - {P2  - V22)X22 

Jo 

+Ai(xii  + a;i2  + Q!^  — Cl)  + A2(2:2i  + ^22  + <^2  ~ ^2)  (3-7) 

Partially  differentiating  the  Lagrangian  function  (3.7)  with  respect  to  all  variables,  the 


necessary  conditions  become 
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(/ii +Pi)$£)i(xji  + a;2i)  - (pi  - wii)  + Aj  = 0 

{hi+Pi)i^DiiXu+X*2^)-{pi-V2l)  + X*2  = 0 

(/l2  +^2)^02(3^12  +^22)  - (P2  - -^12)  + Ai  = 0 

(/l2  + ^2)*^* 02  (3^12  + 3^22)  “ (P2  — '^22)  + A2  = 0 

3;ti  + 3:^2  + - Cl  = 0 

3^21  + 3^22  + ^2^  — C2  = 0 


where  (Xj*^,q;*)  is  a local  minimum. 

The  Lagrange  multipliers,  Ai  and  A2,  can  be  either  zero  or  a positive  value  while 
Ofi  and  0:2  can  be  either  zero  or  nonzero.  It  should  be  noted  that  for  any  i,  at  least  Aj 
or  Oi  has  to  be  zero.  For  each  warehouse,  there  are  three  possible  cases  of  Aj  and  ai 
combination.  Thus,  there  are  nine  cases  for  the  two-warehouse  problem  to  be  considered. 
The  necessary  conditions  for  these  cases  are  summarized  in  Table  3.1.  The  optimal 
solution  satisfies  conditions  in  one  of  these  cases.  Suppose  that  the  optimal  At,  A^  ,al 
and  o;2  are  found,  we  obtain  the  solution  from  solving  the  following: 


2X\al  = 0 


2X*2a*2  = 0 
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Table  3.1:  Summary  of  necessary  conditions  for  an  optimal  solution  to  a two- warehouse 
two-retailer  problem 


Case 

Ai 

ai 

A2 

«2 

Conditions 

1 

0 

0 

0 

0 

Vn  = V2i’,  V12  = V22] 

Xii  -f-  X21  -\-  X\2  + 2:22  = C\  -\-  C2 

2 

0 

0 

+ 

0 

Vn  - V21  = V12  - V22;  vn  > V2UV12  > U22; 
Xu  -)-  2:21  4-  2:12  + 2:22  = Cl  + C2 

3 

0 

0 

0 

-/+ 

^^11  = V21]  V12  = V22', 

2:11  4-2:12  = Cl]  X21  4-  2:22  < C2 

4 

0 

-/+ 

0 

0 

^11  = t^2i;  '^^12  = "^22; 

2:11  4-  2:12  < Ci;  2:21  4-  X22  = C2 

5 

+ 

0 

0 

0 

V21  — Uii  = V22  — V12]  U21  > Uii;  V22  > V12] 

Xu  4-  2:21  4-  2:12  4-  X22  = Cl  4-  C*2 

6 

0 

-/+ 

0 

-/+ 

Vu  = V21;Vi2  = V22]Vii  > V21] 
^11  4-  2:21  < Ci;  2:12  4-  2:22  < C2 

7 

0 

-/+ 

+ 

0 

Vu  ~ V21  — Vi2  — V22]  I'll  > U21;  U12  > U22; 
2:11  4-  2:21  < Ci;  2:12  4-  X22  = C2 

8 

+ 

0 

0 

-/+ 

V2I  - i'll  = V22  - V12]  V2I  > Vii]V22  > V12] 
^11  4-  X21  = Ci;  2:12  4-  2:22  < C2 

9 

+ 

0 

-f 

0 

i’21  ~ i’ll  ’ i'22  ~ i’12; 

2:11  4-  2:21  4-  2:12  4-  X22  — Cl  C2 

^12  T ^22  — ^ 


1 

D2 


fP2-  Vl2  - \ 

V h2  + P2  J 


Xi2  -|-  X22  — *1* 


1 

D2 


f P2  - V22  - A;\ 

V h2+P2  / 


xti  -F  x*2  + af  = Cl 


^21  + ^22  + ^2^  — ^2 


3.5.1  Numerical  Example 

Now  suppose  that  Ri  and  R2  have  negative  exponential  demand  distribution  with 
mean  Di  = 100  units/month  and  D2  = 50  units/month.  The  ratios  of  holding  and 
penalty  cost  at  Ri  and  R2  are  given  as  ^ = 20  and  ^ = 25.  The  holding  cost  is 
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evaluated  hy  h = iV,  where  i is  interest  rate  (cost  of  capital)  and  F is  a value  of 
commodity  per  unit.  Thus,  hi  = iiV  and  /i2  = i2V,  where  ii  = 20  and  22  = 15$/$(?/ear). 
The  transportation  cost  per  unit  is  given  as  follows: 

uii  = 10  $/unit;  V12  = 20  $/unit;  V21  = 5 $/unit;  V22  = 15  $/unit 

Finally,  the  capacity  at  both  warehouses  are  equal  to  100  units.  From  the  given  data 
above,  the  Lagrangian  function  becomes 

F(a:ii,  a;i2,  X2i,  X22,  Ai,  A2,  cti,  0:2)  = 


Al(Xii  + X12  + — 100)  + A2(x21  + X22  + <^2  ~ 100) 

There  exist  Lagrange  multipliers  A^  = 0 and  A2  = 5 such  that 

(69.9) $d,(x*i  Tx^i) -56.6  + A*  = 0 

(69.9) $o,(x*i +x;j -61.6  + A;  = 0 
(39)$D2(2;l2  + a;;2)  - 17.5  + At  = 0 
(39)$D,(xt2  + x-;2) -22.5  + A;  = 0 


(39)  / ^nA^2)d^2  + (37.5)(50)  - (37.5)(xi2  + X22)  + (20)xi2  + (15)x22  + 


(69.9)  / $d,(6)c?6  + (66.6)(100)  - (66.6)(xu  + X21)  + (lO)xn  + (5)x2i  + 


Xii  + Xj2  T cni^  — 100  — 0 


X21  T X22  T 0^2^  — 100  — 0 


2X{al  = 0 


0 
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Then  we  obtain  the  following  solutions: 

= 170 
Xi2  + 3:22  = 30 
^11  + ^12  — 100 
X21  + 2:22  — 100 

Solving  this  system  of  linear  equations,  we  obtain  many  values  of  x^2,  X21,  X22  but 
they  all  give  the  same  objective  value  of  $4, 348.14.  The  example  of  solutions  are  = 
70,  x*2  = 30,  x*2i  = 100,  x*22  = 0),  (x*i  = 80,  x^2  = 20,  = 90,  x*2  = 10),  etc. 

3.6  The  Single- Warehouse  Multi-Retailer  Problem 
This  section  considers  a system  consisting  of  one  warehouse  and  n retailers.  The 
problem  can  be  formulated  as  follows: 

min  {hi  + pi)  ^ (^i)d^i  + Pi-Di  — piXi  + ViXi  + 

{hn  T Pn)  Jq  ^D„  (^n)d^n  T Pn^n  Pn^n  T '^n^n 
S.t.  Xi  + . . . + X„  < I 

Xj  > 0 forj  = 1, . . . ,n 

Note  that  since  we  have  only  one  warehouse,  the  subscript  is  suppressed  to  j (index  on 
retailer).  Converting  the  capacity  inequality  constraints  to  the  equality  case,  we  obtain 


Xi  + . . . + x„  + Q;^  — C 
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The  Lagrangian  function  of  the  problem  can  be  expressed  as 


F{xi,...  ,Xn,X,a)  = 


(hi  + Pi) 


Jo 


+ piDi  - piXi  + viXi  + 


rXn 

(hji  ~b  Pn)  I T Pn^n  Pn^n 

Jo 

A(a;i  + . . . + Xji  + — I) 


+ VnXn  + 


The  necessary  conditions  become 


(hi  + pi)^Di  (a:t)  - (pi  - t^i)  + A*  = 0 


(h„  + Pn)^D„  «)  - {Pn  -Vn)  + X*  = 0 


X^-^  + X-^2  T Oi^  — Cl  — 0 


2A*tt*  = 0 


Then  the  solution  is 


X*  = ) forj  = l,...,n 


hj  + Pj 


X*  = arg  min  < > 4> 

A>o  I ^ 


i fPj-  ^3  - 


V hi  + 


< C 


.j=i  '' 

Table  3.2  summarizes  the  possible  values  for  A and  a and  necessary  conditions  for  an 


optimal  solution. 
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Table  3.2:  Summary  of  necessary  conditions  for  an  optimal  solution  to  a single  warehouse 
n retailers  problem 


Case 

A 

a 

Conditions 

1 

0 

0 

hj+pj^ 

2 

0 

-/+ 

Ej^j  < 1 

3 

+ 

0 

(x*)  = 2lZliI±- 

hj+pj  ’ 

T.i^j  = i 

3.7  The  Multi- Warehouse  Multi-Retailer  Problem 
This  section  analyzes  a general  case  when  the  system  is  consisting  of  m warehouses 
and  n retailers.  The  formulation  is  shown  below. 

n ( m ^ n m 

min  <{hi+  Pi)  ^Dj  {Qd^j  + PjDj  ~ Y f + S 

j=l  I t=l  J j=l  i=l 

n 

S.t.  Y^^ij  ^ 

J=1 

Xij  >0  Vz,j 

Converting  m inequality  constraints  to  the  equality  ones,  the  Lagrangian  function 
can  be  expressed  as  follows: 


^Dj{^j)d^j  + PjDj  'y  ^pjXjj  ^ + 

j=i  1=1 


m n 


+ Yj  ^b  + ~ ^i) 

i=l  j=l 
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Equivalently, 


m 


+ — Ci) 


Thus  the  necessary  conditions  become 


m 


(ft,- - (Pj  - Oi,  - A‘)  = 0 Vi,,- 


(3.8) 


n 


J2x’ij  + af-Ci  = 0 Vi 


(3.9) 


2A*q‘  = 0 V, 


(3.10) 


As  mentioned  earlier,  for  any  i at  least  or  a,  has  to  be  zero.  For  each  warehouse, 
(Aj,  ai)  can  be  one  of  three  cases:  (0,  0),  (0,  +/-),  or  (+,  0).  Therefore  for  m warehouses, 
the  total  sets  of  A^  and  oij  combination  are  S'"!  It  is  increasing  exponentially  with  the 
number  of  warehouses. 

Suppose  that  the  optimal  Lagrange  multipliers,  A*’s  and  o;*’s  are  known,  then  the 
solution  can  be  obtained  from 


However,  obtaining  the  A*  and  a*  is  not  trivial.  Thus,  analytically  solving  the  multi- 
retailer multi-retailer  system  is  very  tedious  work. 


(3.11) 


n 


(3.12) 
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3.8  Conclusions 

This  chapter  develops  the  mathematical  model  to  minimize  the  sum  of  inventory 
(holding  and  penalty  costs)  and  transportation  costs  for  a multi- warehouse  multi-retailer 
system  with  stochastic  demand.  The  analysis  initially  considers  the  linear  transportation 
cost. 

We  first  investigate  a simple  model  with  two  warehouses  and  two  retailers  in  the 
system.  The  Lagrange  multiplier  method  is  employed  to  obtain  the  optimal  solution. 
The  optimal  total  amount  received  at  the  retailer  can  be  derived  from  the  ‘newsboy’  for- 
mula assuming  that  the  optimal  Lagrange  multipliers  are  known.  Obtaining  the  optimal 
multipliers  is  complicated.  Thus,  we  derive  a set  of  feasible  multipliers  and  its  necessary 
conditions  to  ensure  optimality.  Then  we  apply  the  same  technique  to  a single  ware- 
house multi-retailer  system.  Lastly,  we  extend  the  result  to  the  m- warehouse  n-retailer 
system.  For  this  system,  the  possible  values  of  multipliers  are  increasing  exponentially 
with  the  number  of  warehouses  in  the  system.  Thus,  the  Lagrange  multiplier  method  is 
not  efficient  for  solving  large  problems.  In  a later  chapter,  we  develop  a more  powerful 
method  which  can  efficiently  and  optimally  solve  the  linear  transportation  cost  problem. 


CHAPTER  4 

THE  FIXED  CHARGE  TRANSPORTATION  COST  MODEL 

4.1  Introduction 

The  previous  chapter  analyzes  the  case  where  we  have  a linear  transportation  cost. 
However,  this  linear  type  of  cost  structure  is  not  practical.  In  this  chapter  we  consider  a 
more  complex  cost  structure,  a fixed  charge  transportation  cost  function,  which  exploits 
economies  of  scale.  The  fixed  charge  cost  function,  which  has  a discontinuity  at  the 
origin,  is  a special  case  of  concave  cost  structure.  Thus,  the  objective  function  has  a 
convex  (inventory  cost)  and  a concave  component  (fixed  charge  transportation  cost) 
which  makes  the  problem  nonconvex  optimization  problem. 

As  mentioned  earlier,  the  fixed  charge  problem  is  a special  case  of  a concave  mini- 
mization problem.  The  early  research  in  minimum  concave  cost  network  flows  problem 
was  done  by  Zangwill  [60].  One  approach  for  this  class  of  problem  is  an  enumerative 
method  by  ranking  the  extreme  points  [46].  This  method  exploits  the  concave  function 
property  that  every  local  and  global  solution  occurs  at  the  extreme  point  of  the  feasible 
domain  [48,  49].  Implementing  this  approach  to  a large  scale  problem  is  not  practical 
because  of  the  exhaustive  use  of  the  computational  resource. 

The  most  utilized  approach  for  the  fixed  charge  problem  is  branch  and  bound  tech- 
nique. An  exact  solution  is  found  by  transforming  the  fixed  charge  problem  to  a zero-one 
mixed  integer  programming  (MIP).  The  branch  and  bound  method  finds  an  exact  solu- 
tion but  it  takes  a lot  of  computational  efforts. 
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Alternatively,  many  studies  focus  on  developing  heuristics  to  estimate  the  optimal 
solution.  One  of  the  most  efficient  heuristic  is  the  Dynamic  Slope  Scaling  Procedure 
(DSSP)  developed  by  Kim  and  Pardalos  [34,  35].  This  heuristic  is  tested  on  a wide 
range  of  network  structured  problems.  The  reported  computational  results  show  that 
the  DSSP  obtains  solutions  optimally  or  near  optimality  within  reasonable  time.  The 
heuristic  can  approximate  the  optimal  solution  for  problems  as  large  as  10000  arcs  in  a 
little  more  than  one  minute  of  CPU  time. 

The  material  of  this  chapter  is  organized  along  the  following  sequence.  Section  3.2 
approximates  the  inventory  cost  component  using  Monte  Carlo  method,  and  incorpo- 
rates a fixed  charge  transportation  cost.  Section  3.3  discusses  the  DSSP  to  estimate 
the  optimal  of  an  approximated  problem.  Section  3.4  focuses  on  the  computational 
experiment.  The  conclusions  are  presented  in  Section  3.5. 


4.2  Model  Formulation 

4.2.1  Approximated  Problem 

Recall  the  integrated  inventory  and  transportation  problem  formulation. 

n m 

min  {hj  fp=^  CY  Xij  - {Qd^j+ 

i=\ 

m 

Pj  Xij  ~ i^j)d^j)  + 


i=i 


1=1 


'Yj  'Y/ 

j=i  i=i 
n 

s.t.  Y^^n  — 
1=1 


X. 


V*,  j 


> 0 
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Algorithm  IT-1 

' 1.  Generate  U from  Uniform{0,l) 

2.  X ^ F^\U) 

3.  Deliver  X 


Figure  4.1:  Inverse  Transformation  (IT)  algorithm 

The  most  complicated  component  in  the  objective  function  above  is  the  integrand  of 
inventory  cost  component.  In  order  to  eliminate  the  integration  while  trying  to  maintain 
the  randomness  property  of  the  demand  function,  we  would  estimate  the  integral  by 
using  numerical  integration  of  Monte  Carlo  method.  The  Monte  Carlo  method  changes 
the  integration  problem  into  a statistical  estimation  problem.  In  the  changing  process, 
first  we  follow  the  inverse  probability  integral  transformation  method  to  get  sample 
points.  Then  the  integral  is  estimated  by  the  sample  mean.  Figure  4.1  presents  the 
inverse  transform  method,  where  A is  a random  variable  with  cumulative  probability 
distribution  function  (c.d.f)  Fx{x)  [50].  The  algorithm  needs  one  random  number  to 
generate  a sample  point.  Then  this  random  number  is  transformed  to  a sample  point 
through  the  inverse  function.  By  /T  - 1 algorithm,  we  can  randomly  generate  a set  of 
demand  points  rather  than  using  demand  distribution  function.  Then  we  can  estimate 
the  integral  with  the  mean  of  these  generated  demand  points. 

Let  be  a demand  random  sample.  Following  7T  - 1 algorithm,  we  can  obtain 
from  the  equation  below. 


where  $(•)  is  a demand  distribution  function,  and  [/  is  a random  number  from  Uni- 
form(0,l).  The  index  s indicates  the  random  sample  generated  or  we  call  it  scenario 
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for  the  rest  of  this  chapter.  Now,  let  S'  be  a total  number  of  scenarios  generated  to 
estimate  the  integral  of  inventory  cost  component.  We  substitute  the  inventory  cost 
component  with  the  integral  estimating  term,  which  is  the  sample  mean  of  scenarios. 
The  inventory  holding  cost  term 


ho 


rET=i^ij  ^ 

I ( / . ~ ^j)^Dj{^j)d^j 


is  replaced  with 


hi 


S m 


^j) 


s=l  i=l 


Note  that  the  expression  /(•)'^  means  max(0, /(•)).  Thus  if  the  total  of  shipment, 

m 

is  greater  than  the  generated  demand  scenario  (^|)  then  there  occurs  an  excess 
inventory.  In  another  way  around  where  the  generated  demand  scenario  is  greater  then 

TTi  ^ S m 

the  total  shipment,  the  term  would  be  zero.  Thus  — 

i=l  ^ s=l  i=l 

means  an  average  holding  inventory  over  S scenarios  generated.  This  replacing  term  is 
more  attractive  than  the  integration  one.  In  the  same  manner,  the  inventory  shortage 
cost  term 


is  replaced  with 


Pj 


i^j  ^ ^ ^ij )PDj  {^j ) d^j 


1 = 1 


El 

s 


S m 


Then  the  objective  function  becomes 

Ti  S j m n S m n m 

E E 5 E + E E I K]  - E + E E (‘‘■c 

j=l  S=1  i=l  j=l  s=l  i=l  j=l  i=l 
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m m 

There  are  two  complex  terms,  ^ij-CjV  and  {Cj  in  the  expression  above. 

1=1  j=i 

To  simplify,  we  reformulate  the  problem  by  introducing  two  additional  parameters, 
and  Zj~,  and  three  more  constraints: 


> 

iei 

(4.2) 

zr 

> 

(j  — ^ Xij 

iei 

Vi,s 

(4.3) 

> 

0 

(4.4) 

m 


m 


The  terms  and  iCj-Yl  Xij)~^  in  problem  (4.1)  are  replaced  with  and 

i=l  i=l 

Zj~.  These  new  two  variables  represent  the  holding  inventory  and  inventory  shortage 
amount  respectively  which  are  determined  by  the  relation  in  constraints  (4.2)  and  (4.3). 
The  constraint  (4.4)  is  added  to  indicate  that  Zj~^  and  Z^~  always  have  positive  value 
or  zero.  It  is  obviously  to  see  that  Zj~^  and  Zj~  cannot  be  positive  at  the  same  time. 
For  example,  if  the  total  shipment,^  Xij,  is  greater  than  the  scenario  demand,  then 


i€l 


Zj~^  will  have  a positive  value  and  ZJ  is  forced  to  be  zero  by  constraint  (4.4)  through 


the  minimizing  process  in  problem  (4.1). 


Now  the  original  problem  can  be  expressed  as  follows: 


n S , n S n m 

j=l  s=l  j=l  s=l  j—l  i=l 

n 

s.t.  E Ci  Vi 

j=i 

m 

i—\ 

m 


Xij,Z‘+,Z’f  > 0 


Vi,J,s 
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It  should  be  noted  that  the  inventory  cost  component  is  now  a linear  function. 

4.2.2  Fixed  Charge  Transportation  Cost 

In  this  section,  we  will  study  a special  case  of  the  concave  transportation  cost.  As- 
sume that  Tij  is  the  fixed  charge  transportation  cost  function  and  it  can  be  shown  as: 


where  is  a fixed  cost  incurring  whenever  there  is  a flow  from  Wi  to  Rj.  The  Vij,  is  a 
unit  variable  cost  of  transporting  commodity  from  Wi  to  Rj . 

We  can  transform  this  fixed  charge  transportation  problem  as  a zero-one  mixed 
integer  problem  (MIP)  by  introducing  binary  variables  yij.  The  transportation  cost  Tij 
has  a form: 


with 


0 if  Xij  = 0 


1 if  Xij  > 0 


Then  our  problem  becomes: 


mm 


n 


m 


'^Xij  < z^^  + Cj  yj,s 
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16/ 

> 

q-zr 

Vj,s 

Xij 

< 

CiUij 

^ i 

> 

0 

Vij 

— 

{0,1} 

V/,  j 

4.3  Modified  Dynamic  Slope  Scaling  Procedure 
The  Dynamic  Slope  Scaling  Procedure  (DSSP)  was  developed  by  Kim  and  Parda- 
los  [34,  35].  This  heuristic  is  an  efficient  and  effective  solution  approach  for  the  fixed 
charge  and  concave  piecewise  network  flow  problems.  It  can  estimate  the  optimal  so- 
lution with  a small  distance  from  the  optimal  solution  in  reasonable  time.  The  DSSP 
approach  does  not  attempt  to  solve  the  problem  as  a mixed  integer  problem  but  rather 
use  a sequence  of  linear  problems  to  estimate  the  optimal  solution.  The  fixed  charge 
cost  is  represented  as  a “linear”  cost  factor,  where  its  relation  is  expressed  as  follows. 

F 

- t ^ '^3 

— Ujj  “1“ 

Xij 

This  linear  cost  factor  reflects  the  slope  of  the  cost  function  of  linear  problem.  From  the 
relation  above,  Vij  changes  as  Xij  changes.  The  procedure  continues  updating  Vij  until  Xij 
can  make  the  linear  factor  effectively  represents  the  combination  of  variable  cost  Vij  and 
marginal  fixed  cost  Fij  simultaneously.  For  example,  at  point  Xij  in  Figure  4.2,  it  shows 
that  the  solution  of  linear  problem  is  equal  to  the  solution  of  the  fixed  charge  problem 
where  the  heuristic  terminates.  The  computational  test  results  reported  in  [34,  35]  on  the 
problem  size  up  to  37  nodes  335  arcs  show  that  the  DSSP  is  efficient  and  reliable.  The 
worst  percentage  difference  between  the  best  DSSP  solution  and  the  optimal  solution  is 
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Figure  4.2:  Dynamic  slope  scaling  factor  Vij 

0.65%.  Moreover,  the  DSSP  finds  the  solution  in  about  4000  times  faster  than  (measured 
by  CPU  time)  an  exact  algorithm.  The  result  also  shows  that  the  DSSP  can  solve  the 
problem  up  to  202  nodes  10200  arcs  in  a little  over  one  minute  of  CPU  time. 

The  following  DSSP  steps  are  modified  to  test  on  our  approximated  problem. 

Step  0:  Initialize  the  “linear  factor”,  h?  , to  obtain  an  initial  solution,  xF. 

Vij  Vij  Vz,j 


Then,  solve  for  the  initial  solution  using  an  LP. 


Step  1:  Update  the  “linear  factor”,  Vi^^^,  at  iteration  k+1 


V 


Vij  + if  x'^j  >0  at  A:  > 0 


Vij, 


uV . 
O’ 


if  = 0 at  /c  = 0 

•'J 


if  Xij  = 0 at  A;  > 1 
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where  r is  the  index  of  the  most  recent  value  of  the  slope  scaling  factor  when 

Xijr  > 0. 

Step  2:  Solve  the  problem  using  an  LP  with  an  updated  cost  vector. 

Step  3:  Observe  the  solution,  If  (indicate  no  further  improvement) 

then  terminates.  Otherwise  go  to  step  1. 

4.4  Computational  Experiments 

This  section  describes  how  we  measure  the  effectiveness  of  the  DSSP  heuristic  to 
estimate  the  optimal  solution  for  the  approximated  problem.  The  DSSP  solution  is 
compared  with  the  optimal  solution  using  Branch  and  Bound  (B  & B)  method.  Note 
that  we  transform  the  approximated  problem  to  the  mixed  integer  programming  problem 
(MIP)  before  implementing  the  B & B algorithm.  The  two  methods,  DSSP  and  B&B, 
are  coded  in  C with  CPLEX  7.0  (callable  library).  The  computational  experiment  is 
performed  on  an  IBM  SP  PowerPC  Platform  two-processor  node  which  is  332  MHz  with 
512  Mb  of  memory.  The  computational  times  spent  on  both  methods  are  compared  and 
reported.  Finally,  we  can  further  evaluate  the  true  objective  value  from  DSSP  solution 
in  the  original  problem.  This  true  objective  value  is  the  upper  bound  for  the  original 
problem. 

4.4.1  Generating  a Test  Problem 
Parameters  are  set  up  as  follows: 

• Demand 

Demands  at  each  retailer  are  stochastic  and  possess  the  same  distribution  func- 
tion. In  this  study,  we  do  experiments  on  two  distribution  functions:  uniform  and 


55 


exponential.  To  generate  demand  scenarios,  first  we  randomly  generate  the  mean 
of  demand  at  each  retailer  from  the  range  (20,  50). 

mean  = Uniform(20, 50) 

Then  we  generate  a random  number,  rnl,  from  Uniform(0,l).  This  random  number 
is  used  in  the  algorithm  IT-1  (Fig.  4.1)  to  obtain  a scenario  demand  point.  Let  the 
cumulative  distribution  of  demand  (CDF)  at  retailer,  Rj,  be  denoted  by  <F(^j)  = 
where  b is  an  upper  bound  of  the  demand.  Following  the  algorithm  IT-1,  the 
scenario  demand  is  obtained  from 

= rnl  * b 

For  exponential  distribution,  let  the  CDF  be  denoted  by  $(^j)  = 1 — 

We  obtain  the  randomly  generated  scenario  demand  points,  as  follows. 

- —ln{l  — rnl)  * mean 


• Fixed  Cost 

The  fixed  cost  occurs  whenever  there  is  a shipment  from  a warehouse,  VFi,  to  a 
retailer,  Rj.  The  fixed  cost  for  each  link  of  warehouse  and  retailer,  Fij,  is  randomly 
created  from  the  range  (10,  30). 

Fij  = Uniform(10, 30) 

• Variable  Cost 

The  variable  cost  is  a cost  per  unit  of  shipping  commodities  from  Wi  to  Rj.  This 
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cost  occurs  with  the  fixed  cost.  In  practice,  the  variable  transportation  cost  mostly 
depends  on  the  distance  between  the  source  and  the  destination.  Thus  we  generate 
the  variable  cost  based  on  the  distance  between  warehouses  and  retailers.  First  we 
randomly  assign  a coordinate  of  every  warehouse  and  retailer.  Then  the  distance 
between  every  pair  of  Wi  and  Rj  is  calculated.  Finally,  the  variable  cost,  is 
given  from  a range  (1,5)  based  on  the  calculated  distance  above. 

• Penalty  and  Holding  Cost 

The  penalty  and  holding  cost  appear  in  the  inventory  cost  component.  In  this 
study,  the  penalty  and  holding  cost  are  30  and  15  $ per  unit  respectively. 

• Capacity  at  the  Wi 

The  only  set  of  constraints  in  our  problem  is  limited  capacity  at  all  warehouses. 
In  the  computational  experiments,  the  capacity  level,  Ci,  are  the  same  at  every 
warehouse.  The  capacity  level  is  determined  directly  to  the  demand  needed  at  all 
retailers  but  inversely  to  number  of  warehouses  in  each  particular  instance.  For 
example,  in  4 warehouses  5 retailers  problem,  the  capacity  at  the  warehouse  tends 
to  be  more  in  high  demand  instances  than  the  low  demand  instances.  In  addition, 
if  the  number  of  retailers  is  fixed,  the  capacity  at  each  warehouse  tends  to  be  more 
in  the  problem  that  has  less  warehouses.  We  set  the  total  capacity  in  each  instance 
at  70  % of  the  demand  range.  We  do  not  want  to  have  too  high  capacity  because 
the  problem  can  become  an  uncapacitated  problem,  which  is  not  in  our  interest. 
For  the  uniformly  distributed  demands  problem,  the  capacity  at  each  warehouse 
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is  calculated  from  the  equation  below. 


Ci  = ,Dj  + 0.2  * n*  (half  width  of  demand)  | /m 


where  Dj  is  the  mean  of  demand  at  Rj.  The  capacity  level  for  the  exponentially 
distributed  demands  problem  can  be  calculated  from  the  equation  below. 


4.4.2  Computational  Test  Results 

First  we  construct  the  test  problems  in  three  sets  of  various  warehouse-retailer  com- 
binations. The  three  sets  are  two  warehouses  with  various  number  of  retailers,  four 
warehouses  with  various  number  of  retailers,  and  six  warehouses  with  various  number  of 
retailers.  For  each  combination  in  every  set,  we  randomly  generate  10  problem  instances 
for  20,  40,  60,  80,  and  100  demand  scenarios  respectively.  Thus  the  DSSP  heuristic  is 
tested  on  50  problem  instances  for  every  warehouse-retailer  combination. 

The  computational  results  are  reported  in  Table  4.1  - 4.6.  The  first  three  columns, 
W,  R,  S,  indicate  the  number  of  warehouses,  retailers,  and  scenarios  respectively.  Then 
the  following  columns,  we  report  the  average,  minimum,  and  maximum  value  of  relative 
errors  (measuring  the  quality  of  the  solution),  and  CPU  times  in  seconds  using  DSSP 
heuristic  and  B&B  algorithm  respectively.  The  relative  errors  are  calculated  from 


n 


Ci  = (/n(0.3)  * (—  ^ Dj))/m 


where  fossp  is  an  objective  value  obtained  from  DSSP  heuristic,  and  /bub  is  an  optimal 


objective  value  obtained  from  B&B. 
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The  results  show  that  DSSP  finds  almost  optimal  solutions  within  very  small  amount 
of  CPU  times.  For  uniform  demand  case  (see  Tabled. 1,  4.3,  and  4.5),  the  overall  of  the 
relative  errors  average  for  all  test  problems  is  about  1.56%,  and  the  worst  relative  error 
reported  is  8.72%.  For  exponential  demand  case  (see  Table  4.2,  4.4,  and  4. 6), DSSP  gives 
better  estimation.  The  overall  average  of  relative  errors  is  a little  over  1%,  and  the  worst 
relative  error  reported  is  4.51%. 

Figure  4.3  to  4.8  show  the  CPU  times  that  are  correlated  with  number  of  scenar- 
ios, number  of  retailers,  and  the  number  of  warehouses.  The  B&B  time  increases  more 
rapidly  than  the  DSSP  time.  As  a result,  the  gap  between  these  two  times  is  increasing 
exponentially  with  those  two  factors.  It  suggests  that  we  can  trade  the  solution  quality 
with  the  computational  time  by  obtaining  a good  solution  from  the  DSSP  within  rea- 
sonable time  instead  of  obtaining  an  optimal  solution  from  B&B  within  longer  period 
of  computational  time.  It  is  also  noted  that  B&B  times  in  exponential  demand  case  are 
much  less  than  the  uniform  one  (especially  in  problems  of  4 or  more  warehouses) , while 
DSSP  times  are  almost  the  same  in  most  cases  for  both  demand  distributions. 


Table  4.1:  Summary  of  computational  results  for  two  warehouses  with  various  number  of  retailers  problems  with  uniformly 
distributed  demands  
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Figure  4.3:  Average  CPU  times  for  two  warehouses  with  uniformly  distributed  demands 


Table  4.2:  Summary  of  computational  results  for  two  warehouses  with  various  number  of  retailers  problems  with  exponentially 
distributed  demands  
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Figure  4.7:  Average  CPU  times  for  six  warehouses  with  uniformly  distributed  demands 
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Figure  4.8:  Average  CPU  times  for  six  warehouses  with  exponentially  distributed  de- 
mands 
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4.4.3  Evaluating  the  True  Objective  Value 

Our  main  interest  is  the  true  objective  value.  The  feasible  set  of  solution  in  the 
approximated  problem  is  the  same  set  in  the  original  problem  since  they  are  bounded 
by  the  same  set  of  constraints.  Thus  the  DSSP  solution  is  also  a feasible  solution  to 
the  original  problem.  We  can  evaluate  the  true  objective  value  by  borrowing  the  DSSP 
solution  and  substituting  it  in  the  original  objective  function. 


where  lxij>o  is  an  indicator  function.  Its  value  is  one  when  Xij  > 0,  otherwise  zero.  Note 
that  this  is  another  form  of  the  objective  function.  It  is  more  convenient  to  evaluate 
the  objective  value  from  this  form.  From  the  computation  experiments,  we  have  DSSP 
solution  for  uniformly  and  exponentially  distributed  demand.  Let  ^ and 


Substituting  the  demand  distribution  in  Equation  4.5,  the  objective  function  becomes 


^Exi^j)  = 1 - e be  the  CDF  of  uniform  and  exponential  distribution  respectively. 


) T PjDj  PjCy  y ^ij)}  T ^ ^ T ^ij^Xij>o) 


m n m 


for  the  uniform  demand  and 


n 


m 


+ Pj)iY^ 


for  the  exponential  demand. 
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4.5  Conclusions 

This  chapter  approximates  the  integral  components,  inventory  holding  cost  and 
penalty  cost  component,  in  the  objective  function  using  Monte  Carlo  method.  The 
solution  of  this  approximated  problem  with  a fixed  charge  transportation  cost  is  esti- 
mated by  DSSP  heuristic.  The  heuristic  is  tested  with  two  types  of  demand  distribution 
function:  uniform  and  exponential.  Experimenting  with  various  combinations  of  ware- 
houses and  retailers,  the  results  show  that  the  heuristic  finds  solutions  faster  than  the 
B & B method  especially  in  large  problems.  To  measure  the  quality  of  the  heuristic 
solutions,  we  find  a relative  differences  between  the  heuristic  and  optimal  solution.  It 
shows  that  the  heuristic  can  solve  the  approximated  problem  optimally  or  near  optimal- 
ity. It  is  noted  that  the  heuristic  provides  smaller  relative  differences  in  exponentially 
distributed  demands  than  the  uniform  ones.  It  means  that  the  heuristic  works  better 
with  the  exponentially  distributed  demands  case.  Once  we  obtain  the  heuristic  solution 
of  the  approximated  problem,  the  upper  bounds  of  the  original  problem  can  be  evaluated 


from  this  solution. 


CHAPTER  5 

THE  LAGRANGIAN  RELAXATION  MODEL 
5.1  Introduction 

This  chapter  extends  the  study  in  the  fixed  charge  transportation  cost  problem  from 
the  previous  chapter.  Thus  far  the  only  known  fact  about  the  problem  is  the  upper  bound 
while  the  optimal  solution  is  still  unknown.  In  this  case  the  quality  of  upper  bound 
obtained  cannot  be  determined.  As  an  alternative,  the  gap  between  the  upper  bound 
and  lower  bound  can  be  used  as  an  indication  of  the  upper  bound  quality.  Thus  the 
goal  of  this  chapter  is  to  compare  the  upper  bound  obtained  from  the  previous  chapter 
with  the  lower  bound  obtained  in  this  chapter.  In  this  study,  the  lower  bound  is  derived 
by  Lagrangian  relaxation  technique,  which  is  one  of  the  most  successful  techniques  in 
providing  the  lower  bound.  It  is  a general  solving  strategy  [2]  and  easily  applied  to 
wide  variety  of  integer  programming  problems.  In  addition,  the  Lagrangian  relaxation 
provides  much  better  bounds  than  the  linear  programming  (LP)  relaxation. 

This  chapter  is  organized  as  follows.  Section  5.2  presents  a literature  review  of  La- 
grangian relaxation  in  many  diverse  areas.  Section  5.3  describes  the  general  Lagrangian 
relaxation  procedure  and  subgradient  optimization.  Section  5.4  provides  the  Lagrangian 
relaxation  model  for  computing  lower  bounds.  Section  5.5  reports  computational  results. 
Finally  the  conclusions  are  presented  in  Section  5.6. 
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5.2  Literature  Review 

The  purpose  of  this  section  is  not  to  attempt  an  exhaustive  review,  but  rather  to 
survey  the  broad  variety  of  Lagrangian  relaxation  ideas  and  several  possibilities  of  their 
applications. 

The  term  “Lagrangian  Relaxation”  is  coined  by  Geoffrion  [27].  The  key  idea  of 
Lagrangian  relaxation  is  to  remove  complicating  constraints  from  hard  integer  program- 
ming problems.  As  a result,  the  problems  could  “easily”  be  solved  [24,  28,  9j.  In 
many  cases,  the  mathematical  models  can  be  decomposed  to  subproblems  when  the 
Lagrangian  relaxation  method  is  applied.  Typically,  the  subproblem  can  be  efficiently 
solved  by  any  existing  methodology  or  algorithm  [2].  For  example,  the  primary  work 
in  this  area  by  Held  and  Karp  [31]  is  done  in  the  traveling-salesman  problem  where  its 
Lagrangian  problem  is  based  on  minimum  spanning  trees.  In  [24],  it  shows  an  evidence 
that  the  Lagrangian  problem  of  the  generalized  assignment  problem  can  be  dealt  with 
knapsack  solution  procedure. 

Fisher  [25]  presents  the  generic  Lagrangian  relaxation  algorithm  which  consists  of 
three  steps:  (1)  construction  of  branch  and  bound  tree,  (2)  adjustment  of  Lagrange  mul- 
tipliers, and  (3)  solution  of  Lagrangian  problem.  The  first  step  provides  the  best  known 
feasible  value  which  is  passed  to  step  2 and  3.  The  last  two  steps  are  iterated  until  an 
iteration  limit  is  reached  or  a better  bound  is  discovered.  Furthermore,  he  commented 
that  the  relaxed  problem  should  be  easier  to  solve  but  not  too  easy.  It  implies  that  the 
choice  of  relaxed  constraints  also  plays  an  important  role  in  the  success  of  the  relaxation. 
On  top  of  that,  there  is  a question  of  how  to  compute  good  multipliers.  The  answer  to 
this  question  is  that  the  most  utilized  and  general  method  is  the  subgradient  method. 
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Finally,  since  some  constraints  are  eliminated  from  the  problem  formulation,  the  La- 
grangian  subproblem  solution  may  not  be  feasible  to  the  original  problem.  Practically, 
there  is  no  general  rule  of  thumb  to  recover  the  feasible  solution  from  the  Lagrangian 
solution.  It  is  a problem-dependent  procedure. 

As  reviewed  in  Fisher  [24],  literature  oriented  to  the  applications  of  Lagrangian  relax- 
ation during  1970s.  The  Lagrangian  relaxation  was  derived  to  variety  integer  program- 
ming problems  including  the  traveling  salesman  problem,  the  general  integer  problem, 
the  location  problem,  the  general  assignment  problem,  etc.  The  computational  experi- 
ence showed  that  the  Lagrangian  relaxation  provided  extremely  sharp  bounds. 

The  Lagrangian  relaxation  has  become  a heavily  used  tool  in  combinatorial  prob- 
lems. In  the  last  decade,  more  and  more  literature  have  devoted  to  prove  the  power 
of  Lagrangian  relaxation  method  in  many  applications.  An  application  in  the  loca- 
tion problem  is  successfully  done  by  Beasley  [10].  His  computational  results  indicate 
the  robustness  of  Lagrangian  relaxation-based  heuristic  in  four  different  location  prob- 
lems: p-median,  uncapacitated  warehouse  location,  capacitated  warehouse  location  and 
capacitated  warehouse  location  with  single  source  constraints. 

Following  Beasley,  Tragantalernsak  et  al.  [57]  propose  six  Lagrangian-based  heuris- 
tics for  the  two-echelon,  single-source,  capacitated  facility  location  problem.  The  major 
difference  in  these  six  heuristics  is  the  sets  of  relaxed  constraints.  The  computational 
results  have  shown  that  the  relaxed  constraints  affect  the  solution  quality.  This  study 
is  a significant  example  of  deciding  which  constraints  should  be  relaxed.  In  their  second 
paper  [58],  they  develop  a branch  and  bound  algorithm  based  on  the  most  efficient  La- 
grangian heuristic.  The  numerical  results  show  that  the  new  algorithm  requires  smaller 
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branch  and  bound  tree  size  and  CPU  time  than  a standard  LP  based  zero-one  integer 
programming  package. 

The  application  of  Lagrangian  relaxation  in  the  lot-sizing  problem  is  done  by  Afen- 
takis  et  al.  [1].  They  consider  this  problem  in  multistage  assembly  systems.  Their  model 
is  formulated  in  term  of  “echelon  stocks.”  Applied  the  Lagrangian  relaxation  technique, 
the  formulation  is  decomposed  into  single  stage  subproblems.  The  Lagrangian  problem 
provides  sharp  bounds  which  are  consequently  used  to  develop  an  efficient  Branch  and 
Bound  algorithm. 

Balakrishnan  and  Graves  [7]  applied  the  Lagrangian  relaxation  to  develop  a com- 
posite algorithm  in  a multicommodity  network  flow  problem  in  which  the  arc  costs  are 
piecewise  linear  concave  function. 

Recently,  Balakrishnan  and  Geunes  [6]  propose  a Lagrangian-based  heuristic  for  the 
flexible  demand  assignment  problem  (FDA)  which  is  formulated  as  a mixed-integer 
program.  The  heuristic  generates  the  feasible  solution  from  the  decomposed  Lagrangian 
subproblem,  which  can  be  solved  as  a knapsack  problem  with  expandable  items  (KPEI). 
Their  computational  results  show  that  the  upper  bound  obtained  is  tighter  than  the 
linear  programing  relaxation. 


5.3  Lagrangian  Relaxation 

This  section  describes  the  Lagrangian  relaxation  procedure  and  the  detail  in  the 


subgradient  algorithm  employed  to  obtain  the  best  lower  bound. 
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To  demonstrate,  consider  the  integer  problem: 

(IP) 


Zip 

= min  cx 

s.t. 

Ax 

< 

b 

X 

e 

z+ 

The  problem  (IP)  has  one  set  of  explicit  linear  inequality  constraints.  The  decision 
variables  x is  a nonnegative  integer. 

Suppose  that  the  problem  (IP)  is  easier  to  solve  if  the  explicit  constraints  are  re- 
moved or  relaxed.  These  relaxed  constraints  are  brought  into  the  objective  function 
with  associated  penalty  cost  which  is  represented  by  the  Lagrange  multipliers,  A.  Thus 
the  Lagrangian  relaxation  or  Lagrangian  subproblem  is  created  as  follows. 

(LR) 

Zlr  = mincx -h  A(>la;  — 6) 

s.t.  X e 

or  we  can  refer  it  as  the  Lagrangian  function 

L{x,  A)  = min{ca:  -t-  \{Ax  — b)  : x e Z^}. 

As  mentioned  before,  the  Lagrangian  subproblem  can  be  solved  by  exploiting  available 
procedures  or  algorithms  for  a given  set  of  Lagrange  multipliers.  By  weak  duality  the- 
orem [51],  the  optimal  objective  value  obtained  provides  a lower  bound  to  the  original 
problem  (IP).  It  implies  that  Zlh  < Zip.  Since  the  primal  problem  is  a minimization 
problem,  we  wish  to  find  the  greatest  lower  bound.  This  can  be  achieved  through  solving 
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the  Lagrangian  dual  or  the  Lagrangian  multiplier  problem  [2],  This  problem  concerns 
with  choosing  the  multipliers  that  provide  the  sharpest  bound.  It  is  defined  as 

(LD) 


Zld  = msLxL(x,  A) 

A>0 

The  problem  LD  is  a max-min  problem.  We  can  use  an  iterative  approach  to  obtain 
the  best  lower  bound  from  this  relaxation.  That  is,  for  a fixed  value  of  A we  can  solve 
L{x,  A).  Then  to  maximize  Zlr,  we  need  to  find  an  improving  direction  of  A and  take  a 
step  in  that  direction.  Then  we  resolve  the  subproblem  again  with  an  updated  A.  The 
most  common  algorithm  used  is  the  subgradient  optimization{see  details  in  [31]).  This 
algorithm  ensures  that  A converges  to  the  value  that  maximize  Zlr. 

5.3.1  Modified  Subgradient  Optimization  Algorithm 

In  this  subsection,  we  modify  the  algorithm  before  implementing  to  our  problem. 
Let  LED  and  UBD  be  a lower  and  upper  bound  on  the  optimal  value  to  the  original 
problem  respectively.  UBD  is  any  feasible  solution.  Let  t and  dt  denote  an  iteration 
counter  and  the  step  length  parameter  at  iteration  t respectively.  The  description  of  the 
algorithm  is  shown  below. 

Step  1: 

/ 

Set  LBD  = — oo 

t = 0 

< 

do  =2 

Initialize  A°  =0 
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Note  that  LBD  is  determined  by  the  best  Lagrangian  dual  value  and  UBD  is 
determined  by  the  best  feasible  solution  found  so  far.  The  step  scale  {dt)  is  halved 
every  fixed  number  of  consecutively  iterations  where  the  LBD  does  not  improve. 
Step  2:  Solve  minL(x,  A)  and  let  x*  be  the  solution  and  set 

lex 

LBD  = max{LBT>,  L(x*,  A*)} 

Step  3:  Check  for  stopping  criteria 

If  A*  - A*~^  < e for  five  consecutive  iterations  or  t = tmax,  stop.  The  current  dual 
solution  is  chosen  as  the  best  LBD. 

Step  4:  Compute  subgradient 

7i  = - Ci  Vz 

i 

Step  5:  Compute  step  size 


T,=d,(UBD-L(^,X))l\\^f 
Step  6:  Update  the  Lagrange  multipliers 

At+i  ^ ^ 

Step  7:  Update  iteration  counter  (t) 

Note  that  unlike  the  generic  subgradient  algorithm,  we  do  not  attempt  to  find  a feasible 
solution  in  the  algorithm  to  update  UBD.  We  expect  to  obtain  the  maximum  lower 
bound  from  the  subgradient  method.  As  a reminder,  the  purpose  of  our  study  is  to 
measure  the  upper  bound  quality  by  comparing  with  the  best  lower  bound  obtained 
from  the  Lagrangian  problem. 
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5.4  Lasransian  Relaxation  Model 
Consider  the  capacitated  inventory  and  transportation  problem. 

(P) 

” rYlTLi  n m 

min  + pj)  / + p^Dj  - pj{J2  Xij)}  + J2Y1  (5-1) 

J— 1 ^ i=l  j=l  i=l 

subject  to 


n 

Xij  < Ci  Vi 

J=1 

> 0 Vi,j 


(5.2) 

(5.3) 


This  problem  (P)  is  difficult  to  solve  computationally.  However,  the  problem  becomes  an 
uncapacitated  problem  and  could  be  “easily”  solved  if  the  set  of  capacity  constraints  (5.2) 
is  eliminated.  Let  Aj  denote  Lagrangian  multipliers  associated  with  the  set  of  capacity 
constraints  (5.2).  Then,  the  Lagrangian  relaxation  problem  {LR)  can  be  expressed  as 
follows. 

(LR) 

min  Y.{{hj  + pj)  / {Qd^j  + P,  A - PjiY  + 

3 = 1 i^l 

n m m n 

y~!  AAu) + y]  ~ 

j=i  i=i  i=i  j=i 

or 


^ P^YY^  ^ ' ' 771 

+ pj)  / (^j)d^j  + pj  A - ^(Pj  - Ai)xij  + 


m m 


i=l 


i=l 


mm 
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subject  to 


Xij>0  yi,j 


The  Lagrangian  function  is  expressed  as  follows. 


L(x,A)  = 


mm 


Pj^j  "y  \)^ij  + y ^ Tij(xjj) 


m 


1=1 


m 


m 


y ^ ^ 0 


(5.4) 


i=l 


For  given  multipliers,  we  can  solve  the  Lagrangian  problem  (5.4)  optimally.  The 
optimal  solution  of  this  Lagrangian  problem  provides  us  a lower  bound  to  the  original 
problem  (P).  However,  it  does  not  guarantee  that  the  chosen  A will  give  us  a good  lower 
bound.  To  maximize  this  lower  bound,  we  can  solve  to  dual  problem  of  (5.4).  Now,  let 
g(A)  be  the  dual  function  of  minL(x,  A).  Then  the  dual  problem  corresponding  to  the 
problem  (5.4)  is 


As  discussed  eariler,  the  subgradient  method  is  employed  to  obtain  the  best  lower  bound. 

It  is  interesting  that  for  given  A^  for  z = 1, . . . , m,  it  is  noted  that  the  term  Ya=i 
in  function  L(x,  A)  is  a constant  at  this  point.  Then  the  problem  LR  is  decomposed  into 
n uncapacitated  single-retailer  problems  (where  we  have  dropped  the  retailer  index  j). 


(LD) 


maxg(A) 


(5.5) 


st.  A > 0. 
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The  decomposed  function  can  be  expressed  as  follows. 


min 

Xi>0 


+ - K)xi  + 


1=1 


(5.6) 


5.4.1  Fixed  Charge  Transportation  cost 

Now  let  us  incorporate  the  fixed  charge  transportation  cost,  Ti{xi),  to  our  model. 
Ti{xi)  is  defined  as  follows. 


Ti{xi) 


0 if  Xi  = 0 

Fi  + ViXi  if  a:j  > 0 


Substitute  Ti{xi)  in  5.6,  we  obtain 


/ rJ:ZiXi  - ^ \ 

mm\^{h+p)J^  <^Di0d^  + pD-J2iP-^i)^i  + Y.^F,  + ViXi)j  (5.7) 

Since  we  encounter  with  a non-differentiable  function,  the  minimum  cannot  be  found  by 
directly  differentiating.  Solving  the  subproblem  above  is  not  trivial.  However,  we  can 
exploit  the  theorem  below  to  obtain  an  optimal  solution  to  the  Lagrangian  subproblem. 
Theorem  1 Let  there  be  unlimited  capacity  at  the  warehouses.  Then,  if  the  transporta- 
tion cost  function  Ti  are  concave  and  Ti{0)  = 0,  there  exists  an  optimal  solution  in  which 
each  retailer  will  receive  a shipment  from  at  most  one  warehouse. 

Proof-  The  total  transportation  cost  is  expressed  as 

min|^ri(a;,)|  = min  {Ti(a;i)  + T2(x2)  + . . . + r^(a:^)}  (5.8) 

Let  Q denote  the  total  amount  required  at  the  retailer.  Then 


'Y^Xi  = Q 

i=l 
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Now  let 


Xi  = aiQ  Vz 


(5.9) 


where  0 < Oj  < 1,  and  = 1.  Substitute  Equation  5.9  in  Equation  5.8, 

mill  = ^^T^{Ti{aiQ) +T2{a2Q)  + . . . + Tm{ajjiQ)} 


Thus,  the  optimal  solution  is  the  least  cost  among  the  m warehouses.  I 

The  theorem  suggests  that  the  shipment  to  any  retailer  Rj  is  coming  from  only 
one  warehouse,  say  Wi-.,  where  the  minimum  cost  occurs.  Thus,  the  solution  to  each 
uncapacitated  single  retailer  problem  (5.7)  is  to  choose  the  minimum  total  expected  cost 
option  among  m+1  options  (m  options  of  ordering  from  only  one  warehouse  out  of  m 
warehouses  and  one  option  from  ordering  nothing).  If  the  solution  is  to  place  an  order, 
then  the  warehouse  i*j  will  send  the  shipment  to  retailer  j. 

To  determine  the  shipping  amount,  we  can  drop  the  fixed  transportation  cost  from 
(5.7)  first.  Then  we  use  the  classical  newsboy  problem  formula  to  obtain  the  shipment 
size.  For  every  warehouse  i,  the  solution  to  (5.7)  can  be  expressed  as  follows. 


(by  concavity  property)  < min  {aiTi(Q)  + a2T2{Q)  + . . . + a^T^(g)} 


< 


mm{Ti{Q),T2{Q),...  ,o^T^(Q)} 


(5.10) 


The  cost  associated  with  this  decision  is 


(5.11) 
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If  the  decision  is  not  ‘not  ordering”,  its  associated  cost  would  be  pD. 

The  procedure  to  solve  minL(x,  A)  is  summarized  as  follows. 

Step  1:  Find  the  most  economical  option  among  m+1  options  for  every  retailer  and  let 
it  be  denoted  by  TC*. 

Tc;  = 

min  {pjDj, 

+ Pj^j  - {Pj  - + VijX*j  + Fij^  ^ 

Note  that  if  the  decision  is  to  order,  then  we  will  get  the  index  i*  and  the  cost 
associated  with  this  decision  in  the  current  step. 

Step  2:  Sum  up  cost  from  all  retailers  and  subtract  with  the  term  KCi- 

m 

L(x,A)=  ^ TC'  + J]A.Q 

j=l,...,n  i=l 

5.5  Computational  Results 

The  objective  of  the  computation  experiment  is  to  obtain  lower  bounds  in  order  to 
measure  the  upper  bound  quality  obtained  by  employing  DSSP  heurisctic  from  the  pre- 
vious chapter.  To  obtain  the  best  lower  bound,  we  implement  the  modified  subgradient 
optimization  algorithm  discussed  in  previous  section.  The  subgradient  algorithm  is  pro- 
grammed in  C.  The  experiments  have  been  performed  on  the  same  set  of  test  problems 
in  the  previous  chapter. 

We  set  the  subgradient  algorithm  stopping  criteria  parameters  as  follows. 

£ = 10“^ 
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or 


tmax  = 3000 


It  means  that  the  algorithm  would  terminate  if  the  differences  between  Lagrangian 
multipliers  is  less  than  10~^  for  five  consecutive  iterations  or  the  number  of  iterations 
exceeds  3,000.  In  addition,  while  both  stopping  criteria  has  not  been  reached  yet,  if  the 
lower  bound  {LBD)  does  not  improve  for  25  consecutively  iterations,  then  we  will  halve 
the  step  scale  (dt). 

To  determine  the  UBD  in  step  5 (for  computing  the  step  size)  of  the  subgradient 
method,  we  assume  that 


- % 

^ hj  + pj 


Xij  - $ 


)/m 


Then  the  UBD  is  calculated  from 


" rT.T=i  _ JIL  " 

UBD  = + Pj)  / ^Dj  {Qd^j  + pjDj  - PjiY^  Xij)]  + "^{Fij  + VijXij) 

J=1  i-l  j=\  i=l 

For  every  warehouse-retailer  combination  problem,  the  upper  bound  is  compared 
with  two  lower  bounds:  the  initial  and  the  best  lower  bound.  The  initial  lower  bound 
which  is  obtained  by  letting  all  Lagrange  multipliers  be  zero,  A = 0.  Then  the  subgra- 
dient algorithm  is  applied  to  obtain  the  best  lower  bound.  For  each  warehouse-retailer 
combination  problem,  there  are  five  upper  bounds  to  be  compared.  Each  upper  bound 
is  an  estimation  to  the  optimal  value  with  with  different  number  of  scenarios  in  the 


approximated  problem. 
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The  computational  results  are  reported  in  Table  5.1  to  ??.  The  first  three  columns 
represent  the  number  of  warehouses  (W),  retailers  (R),  and  scenarios  (S)  respectively. 
The  following  two  columns  report  the  average,  minimum,  maximum  value  of  the  per- 
centage of  the  relative  gap  %GAP  comparing  the  DSSP  heuristic  solution,  fvssp,  with 
the  best  (maximum)  lower  bound  from  the  Lagrangian  relaxation  model  /lr  which  is 
obtained  from  the  subgradient  method.  The  relative  gap  percentage  is  calculated  from 

%GAP  = (^Lbssp  -Ilr^  ^ 

ILR 

The  results  show  that  for  most  of  the  problems  in  both  demand  distributions,  the 
average  %GAP  behaves  inversely  with  the  number  of  scenarios  in  the  approximated 
problem.  That  means  if  the  number  of  scenarios  increases,  then  the  average  %GAP 
decreases.  In  other  words,  DSSP  can  provide  a better  estimation  to  the  optimal  value 
when  the  distribution  demand  is  approximated  with  more  number  of  scenarios.  The 
results  also  indicate  that  the  average  %GAP  goes  with  the  number  of  warehouses.  Fur- 
thermore, it  appears  in  many  problems  that  the  %GAP  of  the  exponential  demand 
distribution  is  smaller  than  the  uniform  demand  distribution.  It  implies  that  the  the 
DSSP  estimates  the  better  upper  bound  for  the  exponential  distribution.  The  overall 
average  %GAP  reported  for  all  test  problems  are  3.75%  and  2.67%  for  the  uniform  and 
exponential  demand  distribution  respectively.  The  largest  gap  reported  is  15.47%  for 
the  uniformly  distributed  demand  and  9.09%  for  the  exponentially  distributed  demand. 
The  smallest  gap  reported  is  0.32%  for  the  uniformly  distributed  demand  and  0.46%  for 
the  exponentially  distributed  demand. 
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Figure  5.1-  5.4  are  the  line  graphs  showing  the  relationship  between  the  average 
%GAP  and  the  number  of  scenarios.  Each  figure  represents  the  problem  of  a particular 
number  of  warehouses  with  various  number  of  retailers  in  the  problem  except  for  the 
large  problem  graph.  In  the  two,  four,  and  six  warehouses  problems,  each  line  in  the 
graph  represents  the  number  of  retailers  in  the  problem.  In  the  large  problem  graph, 
each  line  represents  a particular  problem  size.  There  are  two  graphs  in  each  figure: 
one  for  the  uniformly,  and  another  one  for  the  exponentially  distributed  demand.  All 
figures  show  that  there  is  less  variability  in  the  average  %GAP  among  the  number  of 
retailers  for  the  exponential  demand.  It  is  easily  to  see  that  the  line  graphs  stay  close 
together  in  the  exponential  demand  graphs  compared  to  the  uniform  demand  graphs. 
That  means  the  number  of  retailers  has  more  effect  in  the  uniform  demand  case  than 


the  exponentially  one. 
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Table  5.1:  Summary  of  % GAP  for  two  warehouses  with  various  number  of  retailers 
problems  with  uniformly  distributed  demands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

2 

5 

20 

5.2705 

2.1947 

8.3672 

40 

2.8661 

0.5144 

5.3046 

60 

2.3473 

0.3237 

7.2278 

80 

2.4486 

0.4749 

3.7548 

100 

2.1089 

0.6396 

5.2497 

2 

10 

20 

4.2828 

2.4421 

9.0450 

40 

2.4131 

1.0785 

5.3973 

60 

2.1173 

0.4425 

3.9709 

80 

1.9033 

1.2114 

2.9769 

100 

1.7756 

1.0815 

3.6845 

2 

15 

20 

4.2896 

2.1811 

6.0484 

40 

2.3952 

0.9198 

3.5147 

60 

2.0099 

1.2727 

3.1396 

80 

1.4892 

0.8115 

2.7431 

100 

1.4213 

0.7224 

1.8799 

2 

20 

20 

4.5427 

2.4838 

6.8975 

40 

2.3394 

1.6414 

3.8105 

60 

1.7766 

1.1182 

3.1927 

80 

1.2272 

0.5507 

1.9643 

100 

1.3011 

0.7467 

1.9137 
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Table  5.2:  Summary  of  % GAP  for  two  warehouses  with  various  number  of  retailers 
problems  with  exponentially  distributed  demands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

2 

5 

20 

3.9758 

2.0715 

6.5749 

40 

2.1036 

0.4584 

3.7732 

60 

1.8275 

0.5179 

3.4639 

80 

1.7294 

1.0000 

3.4693 

100 

1.5203 

0.6977 

2.9473 

2 

10 

20 

4.2632 

2.7749 

7.1038 

40 

2.2597 

1.1128 

3.6709 

60 

1.4194 

0.4608 

2.5744 

80 

1.3361 

0.7674 

1.9985 

100 

1.1641 

0.5850 

1.7388 

2 

15 

20 

3.8097 

1.4686 

5.7992 

40 

2.3610 

1.2213 

3.2544 

60 

1.5932 

1.0595 

2.1216 

80 

1.3221 

0.7468 

2.0919 

100 

1.3383 

0.6574 

1.9184 

2 

20 

20 

4.0781 

1.5664 

6.8344 

40 

2.0801 

1.4231 

3.1403 

60 

1.5680 

1.0088 

1.9857 

80 

1.0642 

0.5010 

1.8217 

100 

1.0405 

0.7005 

1.3780 
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Figure  5.1:  Average  %Gap  (best  LED)  for  two  warehouses 
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Table  5.3:  Summary  of  % GAP  for  four  warehouses  with  various  number  of  retailers 
problems  with  uniformly  distributed  demands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

4 

5 

20 

7.5432 

4.2240 

15.4659 

40 

5.3799 

2.6050 

7.8394 

60 

6.0131 

2.6515 

8.9747 

80 

4.5558 

2.6667 

8.5065 

100 

4.2203 

2.3575 

7.2226 

4 

10 

20 

6.6453 

2.3388 

9.5158 

40 

3.8964 

1.7273 

6.6988 

60 

3.8021 

1.3788 

7.1381 

80 

2.8345 

1.4248 

4.1387 

100 

2.5063 

1.8717 

3.0253 

4 

15 

20 

5.6505 

2.8876 

9.2673 

40 

3.5537 

2.3161 

4.9729 

60 

2.2577 

1.2677 

3.3329 

80 

2.3171 

1.3290 

3.0473 

100 

2.2411 

1.3115 

4.2550 

4 

20 

20 

5.6052 

2.9856 

8.8079 

40 

2.8911 

1.8466 

5.5997 

60 

2.9502 

1.4567 

3.7831 

80 

2.4564 

1.5558 

3.6724 

100 

1.8986 

0.7881 

3.6035 
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Table  5.4:  Summary  of  % GAP  for  four  warehouses  with  various  number  of  retailers 
problems  with  exponentially  distributed  demands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

4 

5 

20 

4.43307 

1.61710 

6.71725 

40 

3.74925 

0.66155 

6.95340 

60 

3.01109 

0.94744 

6.00608 

80 

2.05161 

0.49588 

4.32933 

100 

2.20688 

0.97494 

3.31930 

4 

10 

20 

5.43776 

2.81965 

9.08982 

40 

2.25514 

0.84634 

4.81623 

60 

2.55883 

0.80639 

5.55310 

80 

2.19358 

1.38808 

3.20357 

100 

2.02194 

0.98386 

3.10035 

4 

15 

20 

3.92325 

1.45855 

5.92698 

40 

2.88050 

1.96739 

4.30516 

60 

1.87539 

1.44528 

2.53805 

80 

1.67920 

0.47450 

3.34270 

100 

1.46930 

0.67668 

2.29069 

4 

20 

20 

4.45924 

2.58867 

6.79717 

40 

2.11634 

1.04878 

3.74961 

60 

1.95216 

0.99052 

2.96083 

80 

1.44035 

0.69649 

2.26821 

100 

1.47703 

0.77251 

2.31056 
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Figure  5.2:  Average  %Gap  (best  LED)  for  four  warehouses 
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Table  5.5:  Summary  of  % GAP  for  six  warehouses  with  various  number  of  retailers 
problems  with  uniformly  distributed  demands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

6 

10 

20 

6.35337 

4.79593 

8.15290 

40 

6.54585 

3.92919 

10.06927 

60 

4.39570 

3.02571 

6.91162 

80 

4.56873 

3.03198 

6.57606 

100 

3.61513 

1.73449 

5.79110 

6 

15 

20 

6.01000 

2.97233 

9.65109 

40 

4.74029 

2.33175 

6.94994 

60 

3.59946 

2.50648 

5.28953 

80 

3.70633 

2.35134 

5.04463 

100 

3.15663 

1.94358 

4.38206 

6 

20 

20 

5.50709 

2.40467 

7.68615 

40 

4.49438 

3.08043 

5.26264 

60 

3.54298 

2.15763 

5.90731 

80 

3.03590 

1.63801 

4.15084 

100 

2.93561 

2.04458 

5.02530 

6 

25 

20 

5.25893 

3.09068 

8.08560 

40 

4.31816 

3.40861 

5.87747 

60 

3.27204 

1.95002 

4.30819 

80 

2.41209 

1.58404 

3.16097 

100 

2.49047 

1.49162 

3.85256 
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Table  5.6:  Summary  of  % GAP  for  six  warehouses  with  various  number  of  retailers 
problems  with  exponentially  distributed  demands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

6 

10 

20 

5.1844 

3.0374 

7.3217 

40 

4.4792 

3.0666 

6.3712 

60 

2.8532 

1.0308 

4.5031 

80 

2.6185 

1.2381 

5.3927 

100 

2.5839 

1.0358 

3.3324 

6 

15 

20 

4.5123 

2.9377 

7.0596 

40 

2.7071 

0.9780 

4.5773 

60 

2.5131 

1.5962 

3.8610 

80 

2.2996 

1.7520 

3.2678 

100 

2.0474 

1.4501 

3.1322 

6 

20 

20 

4.0925 

2.2568 

6.5534 

40 

3.0229 

1.5473 

4.3658 

60 

2.2977 

1.5256 

3.5423 

80 

1.7468 

1.2222 

2.3836 

100 

1.9640 

1.2248 

3.1099 

6 

25 

20 

4.7347 

2.5655 

7.3597 

40 

2.8712 

2.1149 

4.0509 

60 

2.2124 

1.4247 

3.1610 

80 

1.6719 

0.8583 

2.2778 

100 

1.5077 

1.2253 

2.0017 
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Figure  5.3;  Average  %Gap  (best  LED)  for  six  warehouses 
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Table  5.7:  Summary  of  % GAP  for  large  problems  with  uniformly  distributed  demands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

10 

50 

20 

5.5157 

5.0315 

6.0444 

40 

3.4126 

2.3851 

4.3306 

60 

2.9608 

2.0042 

3.7984 

80 

2.9218 

2.0405 

3.5223 

100 

2.7513 

2.4068 

3.3224 

10 

100 

20 

5.2196 

4.3332 

6.2107 

40 

3.3914 

2.5652 

3.9007 

60 

2.4356 

1.8034 

2.9868 

80 

2.1334 

1.5379 

2.7518 

100 

1.9974 

1.3347 

2.5310 

20 

50 

20 

7.0553 

5.7918 

8.2267 

40 

5.5698 

4.6696 

7.1511 

60 

4.5981 

3.8933 

5.7933 

80 

4.3412 

3.5291 

5.7252 

100 

3.9715 

3.2088 

4.8972 

25 

50 

20 

7.6140 

6.3478 

8.6434 

40 

6.0406 

5.1851 

6.8680 

60 

4.8569 

4.1901 

5.6154 

80 

5.0438 

4.4815 

5.6444 

100 

4.6679 

3.8926 

5.6179 
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Table  5.8:  Summary  of  % GAP  for  large  problems  with  exponentially  distributed  de- 
mands 


Size 

Scenarios 

%GAP 

W 

R 

Average 

Min. 

Max. 

10 

50 

20 

4.8406 

4.1123 

6.7421 

40 

2.8171 

2.0137 

3.9764 

60 

2.2431 

1.8896 

2.8806 

80 

1.8559 

1.4454 

2.2157 

100 

1.6693 

1.3541 

2.3221 

10 

100 

20 

4.4455 

3.7486 

5.3499 

40 

2.4968 

2.0305 

3.0085 

60 

1.8160 

1.4339 

2.2268 

80 

1.5654 

1.0468 

1.8922 

100 

1.2940 

1.0274 

1.5901 

20 

50 

20 

5.6134 

4.0778 

7.2085 

40 

3.5057 

3.0523 

4.2722 

60 

2.8651 

2.1955 

3.4713 

80 

2.5795 

2.0705 

3.5085 

100 

2.2224 

1.5391 

2.6677 

25 

50 

20 

5.7156 

4.5628 

8.5230 

40 

4.1840 

3.7080 

4.7249 

60 

3.2670 

2.5124 

4.1175 

80 

2.7787 

2.0165 

3.8310 

100 

2.7311 

2.0514 

3.4243 
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Figure  5.4:  Average  %Gap  (best  LBD)  for  large  problems 
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5.6  Conclusion 

This  chapter  measures  the  upper  bound  quality  obtained  from  the  previous  chapter 
using  DSSP  heuristic.  The  criterion  used  is  the  gap  between  upper  bounds  {DSSP  so- 
lution) and  the  best  lower  bounds.  The  lower  bounds  are  derived  from  the  Lagrangian 
relaxation  model.  The  relaxed  problem  can  be  decomposed  into  n single-retailer  un- 
capacitated subproblems.  For  a fixed  value  of  Lagrange  multiplier,  the  solution  to  the 
subproblem  can  be  retrieved  by  using  the  classical  newsboy  problem  formula.  Then  the 
subgradient  algorithm  is  employed  to  solve  the  Lagrangian  dual  problem  to  obtain  the 
best  lower  bound.  The  experiments  have  been  performed  on  various  problem  sizes  with 
two  different  demand  distributions  as  in  the  previous  chapter.  For  both  distributions, 
the  number  of  scenarios  affects  the  DSSP  solution  quality.  The  result  shows  that  if  the 
number  of  scenarios  increases,  then  DSSP  can  provide  a better  estimation  to  the  optimal 
value  {%GAP  is  smaller).  The  overall  average  percentage  gap  of  the  exponentially  dis- 
tributed demand  is  little  better  than  the  uniform  one.  Finally,  the  computational  results 
have  demonstrated  that  the  upper  bounds  obtained  from  the  previous  chapter  can  be 
used  to  estimate  the  optimal  solution  to  the  capacitated  inventory  and  transportation 


problem. 


CHAPTER  6 

THE  LAGRANGIAN  RELAXATION  BASED  HEURISTIC 

6.1  Introduction 

This  chapter  presents  a new  approach  for  a fixed  charge  transportation  model  without 
using  the  scenario  generating  technique  as  in  Chapter  4.  This  new  approach  employs  the 
Lagrangian  relaxation  (LR)  model  to  find  an  optimal  solution  to  the  linear  transporta- 
tion cost  problem.  Unlike  the  scenario-based  DSSP  heuristic,  the  Lagrangian  relaxation 
based  DSSP  heuristic  uses  the  true  demand  distribution  instead  of  an  approximation 
of  the  demand  distribution  (using  scenarios).  After  obtaining  an  optimal  solution  of 
the  linear  transportation  cost  model,  the  linear  cost  factor  is  updated  as  in  the  DSSP 
heuristic. 

The  LR-based  heuristic  obtains  the  solution  by  solving  a series  of  linear  transporta- 
tion cost  problems  (nonlinear)  instead  of  solving  linear  programming  problems  (approx- 
imated problems)  as  in  the  scenario  based  heuristic.  The  LR-based  heuristic  is  better 
estimating  an  optimal  solution  because  the  DSSP  heuristic  is  applied  directly  to  the 
original  problem  not  through  the  approximated  problem. 

The  material  of  this  chapter  is  organized  along  the  following  sequence.  Section  6.2 
reviews  the  linear  transportation  cost  model  and  its  Lagrangian  relaxation  model.  Sec- 
tion 6.3  outlines  the  solution  procedure.  Section  6.4  focuses  on  the  computational  ex- 
periment. Finally,  the  concluding  remark  is  presented  in  Section  6.5. 
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6.2  A Linear  Transportation  Cost  Model 
6.2.1  Model  Formulation 

Let  us  revisit  the  formulation  of  a system  of  m warehouses  n retailers  with  a linear 
transportation  cost. 

n / m 

min  hj  CY  Xij  - {Qd^j 

j=l  V i=l 

m 

'^Pj  /eILi  ~ ^ iCj)dCj 

i=l 

n m 

+ Y 

j=l  i=l 


n 


j=l 

Vi 

O 

Al 

Vi,i 

The  objective  is  to  minimize  the  inventory  cost  and  transportation  cost  subject  to  the  ca- 
pacity and  nonnegativity  constraints.  By  rearranging  terms,  problem  (6.1)  is  equivalent 
to 


n ^ m ^ n m 

min  Y \{hi+Pi)  ^Dj  {Qd^j  + pjDj  - Y Pj^ij  f + 5^ 

I z=l  J j=l  i=l 

n 

s.t.  Y^^ij  ^ C'i  Vi 
1=1 


X 


Vi,j 


> 0 
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6.2.2  A Lagrangian  Relaxation  Model 

Let  Aj,  for  i = 1, . . . , m,  be  a Lagrange  multiplier  corresponding  to  the  capacity 
constraint.  The  Lagrangian  function  is  expressed  as  follows. 


We  can  solve  the  dual  problem  of  Lagrangian  problem  (6.2)  to  yield  the  optimal  Lagrange 
multipliers,  At  for  z = 1, . . . , m.  The  optimal  Lagrange  multipliers  are  distinguished  into 
two  sets  by  their  optimal  values:  A-  > 0 and  A*  =0.  Let  Ilr  denote  the  set  of  warehouses 
which  their  capacity  constraints  are  binding,  and  lut  otherwise.  Thus  Irr  = {i  : A*  > 0} 
and  Ilr  = {*  : A*  = 0}.  We  know  that  if  A*  = 0,  the  constraint  is  redundant.  In 
other  words,  Wi  is  not  used  at  full  capacity.  On  the  other  hand,  if  A*  > 0,  then  the 
constraint  is  binding.  That  means 


m 


(6.2) 


n 


Then  the  original  problem  becomes 


mm 


n m 
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n 


t.  y ^ Xij  — Ci 

Vi  G Ilr 

j=i 

IV 

o 

Vi,J 

The  corresponding  Lagrangian  problem  is 


min  F{xij,  Xi)  = 


" I rY:T=o^ii  _ ^ ^ 

/ y ] Tpj)  / ^ pjDj  — / ,PjXjj  T ^ ^ VjjXjj 

j=i  I -^0  i=l  i=l 


+ Xj(y  ^ Xjj  Ci) 

i^lLR  i=l 


(6.3) 


s.t.  Xij  > 0 


or 


mm 


in  F{xij,Xi) 


" I rYT=o^a  _ ^ 

2^  { {hj  + Pj)  / {^j)d^j  + PjDj  - 22  PjXij 

j=i  I -^0  i=l 


+ y ^i)xij  + y ^ 


Vij  Xij 


> - Y. 

i^^LR 


(6.4) 


S.t.  Xij  > 0 


Now  suppose  that  the  nonnegativity  constraint  is  ignored.  Taking  the  derivatives  of 
Equation  (6.4)  with  respect  to  X{j,  we  obtain 


dF(xij,Xi) 

= [hj  + Pj)$D.  (2^  Xij)  - Pj  + Vij  + ATpe/LH)  j 

d i=o 


If  the  nonnegative  constraint  is  not  binding,  i.e.  x*j  > 0,  then 


m 

{hj  + Pj)^Dj  iY  - Pj  + % + ^d{ieiLR}  = 0 

i=0 
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m < I 

d)  ^ ^ _ Pj  Vjj 


i=0 

In  set  Ilh,  the  warehouses  become  uncapacitated.  That  means  there  will  be  at  most 
one  warehouse  from  this  set  for  which  x*j  > 0 for  any  j.  In  other  words,  if  Rj  needs 
to  be  supplied  by  warehouses  in  set  Ilr,  only  one  warehouse  with  the  minimum  unit 
transportation  cost  will  be  the  supplier  of  Rj. 


6.3  Solution  Approach 

This  section  outlines  the  procedure  of  the  Lagrangian  based  heuristic.  This  procedure 
approximates  a solution  for  the  fixed  charge  transportation  model  by  solving  successive 
a linear  transportation  cost  problems  using  the  Lagrangian  relaxation  method.  Then 
linear  cost  factor  is  recursively  updated.  The  procedure  composes  of  three  main  stages. 
The  first  stage  finds  the  optimal  Lagrange  multipliers  and  its  objective  value  using  the 
subgradient  optimization  of  the  linear  transportation  cost  problem.  The  obtained  multi- 
plier optimal  values  reveal  which  capacity  constraints  are  binding  (or  which  warehouse  is 
used  at  full  capacity).  However,  they  do  not  tell  us  the  primal  optimal  value,  x*y  In  or- 
der to  recover  x*j  from  the  Lagrange  dual  problem  solution,  we  employ  the  perturbation 
idea  to  investigate  which  warehouses  are  supplying  a particular  retailer  in  the  second 
stage.  Finally,  the  linear  factor  is  updated  based  on  the  optimal  solution  obtained  from 
the  second  stage. 

After  all  multiplier  optimal  values  are  known,  we  can  compute  the  following  ratios. 


RAij 


Pj  Vij  AT{i£/} 
d"  Pj 


for  yi,j 
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From  these  ratios,  two  cases  can  happen.  In  the  first  case,  for  any  Rj,  all  corresponding 
RAij's  are  different.  These  different  ratios  yield  different  shipment  sizes.  Among  them, 
there  will  be  only  one  shipment  size  which  gives  the  minimum  total  cost.  That  means 
there  is  only  one  warehouse  supplying  this  particular  Rj.  In  the  second  case,  there  occur 
duplicated  ratios.  The  duplicated  ratios  yield  the  same  shipment  sizes.  It  implies  that 
the  Rj  will  potentially  be  supplied  by  all  corresponding  IFi’s  where  their  ratios  are  equal. 

In  the  perturbation  stage,  we  aim  to  find  what  warehouse  is  supplying  a particular 
Rj.  Let  Xij  be  a new  decision  variable,  and  Xij  = Xij  + e,  where  e is  a constant  number. 
Assume  that  Rj  receives  the  shipment  size  of  e from  W-  in  which  can  be  any  warehouse. 
If  we  want  to  check  whether  or  not  Wi  is  supplying  a particular  Rj,  we  would  perturb 
the  Lagrangian  problem  by  appending  a constraint,  X{j  > e.  If  Xij  is  indeed  an  optimal 
flow  {xij  > 0),  then  the  Lagrangian  objective  value  does  not  change.  On  the  other  hand, 
the  objective  value  increases,  if  Xij  is  not  in  the  optimal  solution  {xij  = 0).  This  causes 
by  initiating  an  unnecessary  e flow  in  the  solution.  The  perturbation  of  all  possible  flows 
is  performed  one  flow  at  a time.  Let  i be  a warehouse  where  xu  = e.  The  total  cost 
function  of  each  perturbed  problem  for  any  Rj  is  as  follows: 


[hj  Pj)  j ^Dj{^j)d^j  PjDj  Pji^ij  + X'-)  + {vij  + Xi)xij  + + ^i)^ij 


or 


{xij  + e)  ^ 


Pj  Vij  Aj 

hj  + Pj 


Then 
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— 1 f Pj  ^ij 


It  should  be  noted  that  if  pj  — vij  — Aj  < 0,  then  Xij  = 0 and  the  total  cost  becomes 


PjDj  + (%  + A<)e 


For  the  case  that  all  RAi/s  are  different,  if  the  perturbed  problem  objective  value  is 
equal  to  the  subgradient  objective  value  in  the  first  stage,  it  means  that  x*j  > 0 and  x*^ 
satisfies 

a.  \ — Pj  ~ ~ ^i^{ieiLR} 

As  mentioned  earlier,  retailers  with  duplicated  ratios  will  potentially  have  multiple 
shipments  from  warehouses  with  the  same  ratios.  At  most  one  of  these  warehouses  is  not 
used  at  full  capacity  and  the  rest  is  in  the  set  Ilr-  In  addition,  some  of  these  warehouses 
may  supply  other  retailers  which  have  split  shipment.  That  means  it  is  not  necessary 
that  there  is  only  one  split  shipment  retailer  for  all  warehouses  in  set  Ilr.  However, 
at  least  one  of  warehouses  in  this  set  has  only  one  split  shipment  retailer.  Since  all 
other  shipments  (flows)  are  known,  the  remaining  one  with  the  split  shipment  retailer  is 
calculated  through  the  capacity  constraint.  To  demonstrate,  let  x*-.,  i G Ilr  be  the  split 
shipment  at  any  Wi.  The  remaining  flow  at  that  warehouse  is  calculated  as  follows: 
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Let  xL,  i G Ilr  be  the  remaining  one  flow  from  the  noncapacited  warehouse  of  Rj. 
This  flow  is  given  as  follows: 


Step  1 Solve  a linear  transportation  cost  using  the  Lagrangian  relaxation  method.  In 
this  step,  the  optimal  Lagrange  multipliers  and  the  objective  value  are  determined 
using  the  subgradient  optimization  technique. 

Step  2 Compute  the  ratios,  RAij  for  Vz,j  to  divide  retailers  into  two  sets:  different- 
ratio  retailer  and  duplicated-ratio  retailer. 

Step  3 Perturb  the  Lagrangian  problem  with  a new  constraint 


Then  compare  the  perturbed  problem  objective  value  with  the  subgradient  objec- 
tive value  obtained  from  Step  1.  If  they  are  equal,  then  the  perturbing  flow  is 
greater  than  zero  (i.e.,  Wi  is  in  fact  supplying  Rj).  Thus,  Xtj  > 0 and  Xij  = Xij  + e. 
Otherwise,  Xij  = 0. 

Repeat  Step  3 for  all  possible  flows  {xij),  or  until  the  supplying  warehouse  for  all 


The  Lagrangian  based  heuristic  is  summarized  as  follows: 


Step  0 Initialize  the  “linear  factor”,  u?-. 


retailers  is  found. 
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Step  4 Recover  the  primal  solution,  x*^,  calculate  the  (fixed  charge)  objective  value, 
and  keep  the  best  value  found. 

Step  5 Update  the  “linear  factor”,  at  iteration  k+1 


Vij  + if  x^j  >0  at  A:  > 0 


= 


Viji 


if  x%  = 0 at  A:  = 0 


if  Xii  = 0 at  A:  > 1 

tj  — 


where  r is  the  index  of  the  most  recent  value  of  the  slope  scaling  factor  when 


XijT  > 0. 

Step  6 Observe  the  solution,  x^^^ . If  x\^^  = (indicate  no  further  improvement)  or  if 
the  (fixed  charge)  objective  value  does  not  improve  for  five  consecutive  iterations, 
then  terminates.  Otherwise  go  to  step  1. 


6.4  Computational  Experiments 

This  section  compares  the  effectiveness  of  the  Lagrangian  based  heuristic  with  the 
scenario  based  heuristic  presented  in  Chapter  4.  Two  key  factors  in  comparison  are 
%GAP,  which  is  the  distance  of  the  heuristic  solution  from  the  best  lower  bound,  and 
the  computational  times.  The  Lagrangian  heuristic  is  coded  in  C. 

The  proposed  heuristic  is  tested  on  three  sets  of  warehouse-retailer  combinations: 
two  warehouses,  four  warehouses,  and  six  warehouses  with  various  number  of  retailers. 
For  every  warehouse-retailer  combination,  50  problem  instances  with  uniform  demand 
distribution  are  tested.  These  test  problems  are  generated  based  on  the  parameters 
described  in  Chapter  4.  The  number  of  scenarios  in  the  test  problems  is  100. 


no 


The  computational  result  is  reported  in  Table  6.1.  The  first  two  columns  are  the 
number  of  warehouses  (W)  and  retailers  (R)  respectively.  The  following  six  columns 
report  the  average,  minimum,  and  maximum  value  of  the  percentage  differences  from 
the  best  lower  bound  for  the  Lagrangian  based  heuristic  and  scenario  based  heuristic 
respectively.  Let  /*  be  an  objective  value  attained  from  the  heuristic.  Let  BLB  denote 
the  best  lower  bound  obtaining  from  the  Lagrangian  relaxation  method  described  in 
Chapter  5.  The  percentage  differences  are  calculated  as  follows: 

The  last  six  columns  report  the  average,  minimum,  and  maximum  CPU  times  used  by 
each  heuristic. 

The  results  show  that  the  Lagrangian  based  heuristic  found  almost  optimal  and 
optimal  solution.  It  provides  better  solutions  than  the  scenario  based  heuristic  in  less 
CPU  times  (see  Figure  6.1  and  6.2).  As  the  number  of  retailers  increases,  the  CPU  time 
differences  increases  as  well.  Especially  in  the  two  warehouse  problem,  the  Lagrangian 
based  heuristic  performed  much  superior  than  the  scenario  one.  While  the  worst  %GAP 
reported  is  6.37%  by  using  the  Lagrangian  based  heuristic,  it  is  11.75%  by  using  the 


scenario  based  heuristic. 
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Two  Warehouses 


5 10  15  20 


♦ LR-based 
—o— Scenario-based 


Number  of  Retailers 
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■ ♦ LR-based 
—D— Scenario-based 


Six  Warehouses 
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—o— Scenario-based 


Figure  6.1:  Comparison  of  %GAP  from  the  best  lower  bound 
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Two  Warehouses 


Four  Warehouses 


— LR-based 
—Cl— Scenario-based 


Six  Warehouses 


♦ LR-based 
— □ — Scenario-based 


Figure  6.2:  Comparison  of  CPU  times 
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6.5  Concluding  Remarks 

This  chapter  develops  a new  Lagrangian  based  DSSP  heuristic.  This  new  heuristic 
eliminates  the  scenario  generating  scheme.  It  approximates  the  optimal  solution  from 
the  true  demand  distribution  not  from  an  approximation  of  demand  distribution  as  in 
the  scenario  based  DSSP  henristic.  Based  on  a set  of  test  problems,  the  result  shows 
that  the  new  heuristic  generates  better  solutions  in  less  CPU  time. 


CHAPTER  7 
FUTURE  RESEARCH 

This  dissertation  investigates  the  issues  of  optimizing  an  integrated  inventory  and 
transportation  decision  for  a multi-warehouse  multi-retailer  distribution  system  in  the 
face  of  uncertainties.  The  primary  contributions  of  this  dissertation  are 

• The  development  of  an  optimization  model  that  extends  the  single-period  stochas- 
tic inventory  problem  to  accommodate  a transportation  cost  and  capacity  restric- 
tions. 

• The  investigation  of  the  heuristic  performance  through  the  computational  experi- 
ments. 

The  current  optimization  model  is  developed  for  a single-commodity  single-period 
problem  concerning  only  a distribution  system  in  the  supply  chain.  This  allows  for 
several  model  extensions.  Three  main  directions  for  future  research  are 

1.  Generalize  transportation  cost  structures 

2.  Extend  to  handle  multiple  commodities  and  multiple  time  periods 

3.  Incorporate  the  production  cost 

In  this  dissertation,  only  linear  and  fixed  charge  cost  structures  are  examined.  The 
more  complex  cost  structure  such  as  the  piecewise  linear  concave  has  not  been  addressed. 
This  cost  structure  is  often  found  in  practice  when  all  units  quantity  discount  structures 
are  applied.  The  piecewise  linear  concave  cost  structures  contains  breakpoints  at  which 
the  slope  of  the  variable  cost  per  unit  changes,  denoted  by  bk^j’s  in  Figure  7.1.  The  fixed 
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Figure  7.1:  Piecewise  linear  concave  cost  structure 

charge  cost  structure  is  the  piecewise  linear  concave  cost  with  only  one  slope  interval. 
Thus  the  cost  in  each  interval  can  be  represented  by  a fixed  charge  cost  function.  In 
addition,  an  important  property  of  the  piecewise  linear  concave  cost  is  that  the  unit 
cost  slope  decreases  as  the  interval  contains  a larger  shipment.  For  example,  the  slope 
in  interval  (0,  bk}^)  is  greater  than  the  slope  in  the  interval  {bkk,  bkfj)  and  so  on. 

As  mentioned  before,  the  optimization  model  in  this  dissertation  concerns  only 
single-commodity  single-period  problem.  This  model  should  extend  to  handle  multi- 
ple commodities  and  multiple  time  periods.  Still,  various  issues  within  this  integrated 
inventory-transportation  model  are  worth  exploring.  These  include  positive  leadtimes, 
backlogging,  and  so  on. 

Finally,  the  management  of  the  supply  chain  operations  is  not  complete  without 
the  production  activity.  The  future  research  should  integrate  the  production  planning, 
inventory,  and  transportation  decisions. 
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